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The  major  objective  of  this  research  was  to  develop  new  test  procedures 
to  investigate  the  fracture  toughness  of  laminated  stitched  composites.  Graphite- 
epoxy  laminated  composites  have  very  high  stiff ness-to-weight  and  strength-to- 
weight  ratios  that  make  them  very  attractive  in  structural  applications.  The 
properties  of  these  materials  depend  strongly  on  fiber  orientation  and  ratio  of 
fiber  to  resin.  The  strength  along  the  translaminar  direction  of  this  material  is 
dominated  by  the  matrix  and  is  typically  about  5%  of  the  strength  along  the  fiber 
direction.  Hence,  these  materials  are  very  easy  to  delaminate.  Through-the- 
thickness  stitching  is  one  of  the  ideal  ways  to  improve  the  translaminar  properties 
and  prevent  crack  propagation.  In  order  to  understand  the  stitch’s  effect,  a series 
of  test  procedures  was  developed  and  was  successfully  used  to  test  Mode  I, 
Mode  II  and  Mixed  Mode  fracture  toughness  of  stitched  laminated  composites. 
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A secondary  objective  was  to  model  stitched  laminated  composites  to 
predict  their  properties.  Finite  element  models  were  proposed  to  understand  the 
effect  of  stitching  on  improving  the  fracture  toughness  for  Mode  I and  Mode  II  of 
stitched  laminated  composites.  From  these  models,  “effective  G”  (effective  strain 
energy  release  rate)  and  “parent  G”  (strain  energy  release  rate  of  the  parent 
laminate)  can  be  isolated  where  stitches  can  be  treated  as  “an  extra  structure.” 
The  “parent  G”  is  an  intrinsic  parent  material  property  and  almost  never  changes 
regardless  of  whether  the  laminated  composites  are  stitched  or  not.  However, 
the  “effective  G”  is  dramatically  increased  by  the  stitching.  The  author  also  found 
that  changing  the  thickness  of  stitches  has  a more  profound  impact  on  the 
“effective  G”  than  changing  the  density  of  stitches. 
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CHAPTER  1 

BACKGROUND  AND  LITERATURE  SURVEY 
1.1  Background 

Graphite/epoxy  composite  laminates  have  very  high  strength-to-weight 
and  stiffness-to-weight  ratios.  They  have  tremendous  advantages  over 
conventional  materials  in  applications  such  as  aerospace  and  automotive 
structures,  in  which  these  two  parameters  are  critical.  Also,  composite  materials 
can  easily  be  tailored  to  obtain  desired  properties  in  different  directions  and  can 
be  optimized  to  meet  specific  performance  requirements.  There  are  other 
advantages  such  as  superior  corrosion  resistance,  high  energy  absorption,  low 
coefficient  of  thermal  expansion,  favorable  thermal  insulation  and  electrical 
resistivity.  These  materials  have  the  potential  to  revolutionize  the  way  present 
day  structures  are  made. 

However,  despite  the  numerous  potential  applications,  production  volumes 
have  not  increased  significantly  in  recent  years.  The  primary  deficiencies  of 
laminated  composites  are  their  poor  interlaminar  strength  and  fracture 
toughness,  low-impact  damage  resistance  and  impact-damage  tolerance,  and 
high  processing  cost.  For  example,  the  tensile  strength  of  carbon  fiber-reinforced 
polymer  composites  is  between  about  500  and  800  MPa  in  the  fiber  direction;  but 
in  the  translaminar  direction  it  is  only  20-30  MPa  since  the  load  is  predominately 
carried  by  the  resin  matrix.  Because  of  the  lack  of  translaminar  reinforcement, 
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they  generally  are  vulnerable  to  delamination  which  is  an  interface  crack  or  a 
debond  between  two  adjacent  layers.  Delamination  can  be  initiated  during 
imperfect  manufacture  or  simply  caused  by  the  impact  of  a foreign  body  during 
service.  Since  the  crack  is  inside  the  structure,  it  is  hard  to  detect  by  visual 
inspection.  The  initial  crack  can  grow  rapidly  under  the  subsequent  application  of 
external  load.  The  overall  stiffness  and  loading  capability  of  the  structure  will 
decrease  significantly.  It  may  cause  catastrophic  failure  of  the  composite 
structure.  The  poor  interlaminar  strength  and  fracture  toughness  of  laminated 
composites  have  restricted  their  application  in  many  lightweight  structures. 

Recently,  due  to  adaptation  of  many  textile  technologies,  the  cost  has 
been  dramatically  reduced.  Researchers  have  mainly  focused  on  how  to 
improve  the  interlaminar  properties  of  composites.  Various  techniques  have 
been  considered  to  enhance  the  out-of-plane  properties.  One  of  the  most 
effective  ways  to  increase  interlaminar  fracture  toughness  is  incorporation  of 
translaminar  reinforcements.  Translaminar  reinforcement  (TLR)  has  two  general 
forms:  continuous  and  discontinuous  (Figure  1-1).  Continuous  weaving,  knitting, 
braiding,  rovings,  threads,  yarns  and  tows  can  be  inserted  into  the  lamina  with 
the  use  of  industrial  sewing/stitching  technology.  Discontinuous  translaminar 
reinforcement  (in  the  form  of  short  fibers,  whiskers,  pins,  etc.)  can  also  be  used 
to  bridge  the  interlaminar  region.  The  processes  of  3-D  weaving,  3-D  knitting 
and  3-D  braiding  generally  have  significant  high-volume  fractions  of  fiber  in  all 
three  directions.  These  processes  have  considerably  improved  fracture 
toughness  and  impact  properties,  but  these  methods  reduce  the  proportion  of 
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fibers  along  the  in-plane  direction  and  create  large  resin  pockets  throughout  the 
structure  which  tend  to  deteriorate  the  in-plane  properties.  The  stitching  process, 
on  the  other  hand,  is  unique  in  that  the  stitched  preform  is  not  an  integral 
structure  because  the  through-the-thickness  yarn  is  inserted  into  a traditional  2-D 
preform  as  a secondary  process  after  lay-up.  So,  this  process  avoids  tracking 
large  resin  pockets  inside  the  composite.  There  are  at  least  three  main  reasons 
for  using  stitching:  i)  Cost-effective  method  for  joining  fabric  cloths  by  stitching 
along  the  preform  edges  which  improves  the  ease  of  handling  prior  to  liquid 
molding,  ii)  Effective  method  for  joining  composite  structures  to  provide  high  in- 
plane and  through-thickness  strength  and  iii)  Improve  the  interlaminar  strength 
and  toughness,  and  impact  damage  tolerance.  Low-density  stitching  (number  of 
stitches  per  square  inch  is  less  than  16)  and  high-density  stitching  (number  of 
stitches  per  square  inch  is  greater  than  32)  both  show  significant  improvement  of 
Mode  I fracture  toughness  and  moderate  enhancement  of  Mode  II  fracture 
toughness  while  the  in-plane  properties  only  suffer  a slight  drop. 

Laminated  composites  can  be  stitched  as  either  prepreg  or  preform. 
Unfortunately,  considerable  fiber  damage  is  sustained  during  stitching  prepregs. 
This  tends  to  deteriorate  the  properties  of  the  composite.  On  the  other  hand,  in 
the  resin-free  fabrics  (preform)  a needle  can  be  easily  pulled  through  the  lamina 
with  significantly  reduced  fiber  damage.  The  type  of  industrial  sewing  machine, 
stitch,  needle  and  thread  are  all  factors  that  affect  this  stitching  technique.  There 
are  two  basic  stitching  types  (Figure  1-2).  One  type  of  stitch  (known  as  the  lock 
stitch)  is  mainly  used  in  the  apparel  industry.  It  consists  of  a two-thread  loop 


4 


between  the  needle  and  the  bobbin  threads,  and  requires  access  to  both  the 
bottom  and  the  top  of  the  laminate.  In  the  apparel  industry,  the  intersection  of 
the  bobbin  and  needle  threads  is  concealed  in  the  fabric  structure  for  aesthetic 
purposes.  However  the  thread  intersection  in  the  middle  of  the  fabric  may  cause 
a high  stress  concentration  at  the  intersection  point.  A modified  lock  stitch 
achieved  by  either  tensioning  the  bobbin  thread  or  easing  the  tension  in  the 
needle  thread  may  overcome  this  problem  by  forcing  the  needle  thread  to  travel 
on  the  surface  of  the  laminate.  The  modified  lock  stitch  is  more  damage-tolerant 
since  the  yarn  travels  through  the  thickness  of  the  structure.  The  other  type  of 
stitch  is  the  chain  stitch  which  has  a mechanism  similar  to  that  of  the  lock  stitch. 

1.2  Literature  Survey 

Mignery  et  al.  [1]  investigated  the  use  of  stitching  with  Kevlar  yarn  to 
suppress  delamination  in  graphite/epoxy  laminates.  Results  showed  that  stitches 
effectively  arrested  delamination.  Dexter  and  Funk  [2]  characterized  the  impact 
resistance  and  interlaminar  fracture  toughness  of  quasi-isotropic  graphite-epoxy 
laminates  made  of  unidirectional  Thornel  300-6K  fibers/Hercules  3501-6  resin 
and  stitched  with  polyester  and  Kevlar  yarns.  They  experimented  with  stitch 
parameters  and  found  a significant  drop  in  damaged  areas  of  stitched  laminates 
compared  to  unstitched  laminates  for  the  same  impact  energy.  The  Mode  I 
fracture  toughness,  characterized  by  the  critical  strain  energy  release  rate,  G|C, 
was  found  to  be  about  30  times  higher  for  the  stitched  laminates.  Ogo  [3] 
investigated  the  effect  of  through-the-thickness  stitching  of  plain  woven 
graphite/epoxy  laminates  with  Kevlar  yarn.  The  study  showed  a manifold 
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increase  in  Gic  values  at  the  expense  of  a slight  drop  in  the  in-plane  properties. 
Pelstring  and  Madan  [4]  developed  semi-empirical  formulae  relating  damage 
tolerance  of  a composite  laminate  to  stitching  parameters.  Mode  I critical  strain 
energy  release  rate  was  found  to  be  15  times  greater  than  in  unstitched 
laminates.  Critical  strain  energy  release  rates  decreased  exponentially  with 
increased  stitch  spacing.  Dransfield  et  al.  [5]  investigated  the  effect  of  through- 
the-thickness  stitching  of  plain  woven  graphite/epoxy  laminates  with  Kevlar  yarn. 
Their  study  showed  a manifold  increase  in  Gic  values  at  the  expense  of  a slight 
drop  in  the  in-plane  properties.  Sharma  and  Sankar  [6]  conducted  a study  on  the 
effects  of  stitching  on  interlaminar  fracture  toughness  of  uniweave  textile 
composites.  They  used  the  double  cantilever  beam  setup  to  test  the  low  density 
stitched  specimens  and  the  University  of  Florida  Compression-After-Impact  (UF- 
CAI)  test  fixture  to  investigate  effects  of  stitching  on  sublaminate  buckling 
behavior.  Their  study  showed  that  stitching  has  a profound  effect  on  Mode  I 
fracture  toughness  and  compression  after  impact  (CAI)  strength.  They  also 
found  that  stitching  does  not  increase  the  impact  load  at  which  delaminations 
begin  to  propagate,  but  greatly  reduces  the  extent  of  delamination  growth  at  the 
end  of  the  impact  event.  Jain  et  al.  [7-9]  used  a double-cantilever-beam  (DCB) 
test  for  stitched  composites  and  found  that  the  cantilever  arms  failed  in  bending 
before  crack  propagation.  They  used  aluminum  bars  to  bond  the  surfaces  of  the 
specimen  to  reinforce  the  specimen  to  prevent  the  specimen  arms  from  failing 
during  the  DCB  tests.  A comparison  of  the  specimens  with  tabs  and  without  tabs 
showed  that  tabs  altered  the  failure  mechanism  of  the  stitched  thread;  and  that 
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the  improvement  in  Mode  I delamination  toughness  was  underestimated.  The 
fracture  toughness  of  stitched  composites  had  10-fold  improvement.  Numerous 
studies  [10-20]  used  the  DCB  setup  to  measure  the  Mode  I fracture  toughness  of 
various  stitched  composites  and  found  that  the  DCB  setup  was  effective  for 
unstitched  laminates  or  lightly  stitched  composites  where  the  stitching  density  is 
low  or  the  stitches  are  not  strong.  However,  as  the  stitching  density  increases  or 
the  strength  of  a stitch  increases,  the  DCB  setup  doesn’t  work.  Some 
researchers  tried  to  reinforce  a specimen  with  aluminum  tabs  to  improve  the  load 
capability  of  a specimen.  They  faced  the  same  problems  that  Jain  et  al.  [9] 
encountered  and  found  that  the  tabs  alter  the  basic  mechanism  of  a specimen 
and  the  interfaces  between  the  tabs  and  the  specimen  are  so  weak  that  they 
tend  to  debond  under  high  loading.  After  carefully  studying  the  failure 
mechanism  of  the  stitched  specimens,  Chen,  Ifju  & Sankar  [21-22]  found  a 
conceptual  solution  for  Mode  I testing  of  heavily  stitched  specimens.  Using  this 
idea,  the  authors  developed  a preliminary  fixture  to  successfully  test  heavily 
stitched  specimens.  Chapter  2 gives  more  detail  about  this  process.  Chapter  3 
shows  the  final  design  of  our  loading  fixture  and  some  test  results  and  analysis. 

A study  by  Sharma  and  Sankar  [6]  showed  that  Mode  II  fracture 
toughness  of  stitched  specimens  increases  with  the  delamination  area  bridged  by 
stitches.  They  showed  that  in  Mode  II  loading  stitches  do  not  break  as  widely 
expected,  but  deform  the  surrounding  matrix.  This  deformation  absorbs  some 
energy  from  the  system  and  increases  the  Mode  II  fracture  toughness.  Mode  I 
and  Mode  II  fracture  studies  by  Sharma  and  Sankar  [6]  were  concerned  only  with 
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unidirectional  uniweaves  where  stitching  was  parallel  to  the  fiber  direction.  Their 
study  indicates  that  the  interactions  among  stitches,  fibers  and  matrix  phase  of 
the  parent  laminate  are  quite  complicated.  The  mechanism  has  to  be  thoroughly 
understood  for  proper  selection  of  stitch  materials  and  stitch  density.  Jain  and 
Mai  [25]  used  two  specimen  geometries:  End-Notch  Flexure  (ENF)  geometry  and 
End-notch  cantilever  geometry  to  analyze  the  effect  of  stitching  on  Mode  II 
delaminaiton  toughness.  They  found  that  the  effect  of  stitching  on  both  specimen 
geometries  is  the  same.  Massabo  and  Cox  [26]  used  the  crack  bridging  concept 
to  study  the  effect  of  stitching  in  Mode  II  delamination.  Jain  et  al.  [27] 
investigated  the  Mode  II  delamination  toughness  by  using  an  end-notched  flexure 
specimen  and  found  that  the  application  of  the  load  to  the  unstitched  specimens 
caused  the  initial  delamination  to  grow  in  an  unstable  manner.  With  the  addition 
of  stitching  a similar  loading  situation  resulted  in  stable  crack  growth  from  the 
initial  delamination.  The  improvement  due  to  stitching  was  in  the  range  of  100% 
to  200%  with  respect  to  unstitched  specimens.  It  was  also  observed  that 
increasing  stitch  thread  diameter  or  stitch  density  generally  led  to  larger 
improvement  in  Mode  II  toughness.  Sankar  and  Sharma  [6]  used  three 
approaches  (area  method,  equivalent  area  method  and  classical  beam  theory)  to 
calculate  Gnc.  They  determined  the  location  of  the  crack  tip  with  a microscope 
and  an  ultrasonic  C-scan.  Results  show  that  an  increase  in  apparent  Gnc  of 
stitched  composites  was  5 to  15  times  when  the  crack  was  allowed  to  propagate 
up  to  about  the  center  line  of  the  laminates.  The  stitches  did  not  break  during 
these  tests.  Clara  and  Barry  [28]  investigated  the  accuracy  of  the  four-point 
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bending  end-notched  flexure  (4ENF)  test  for  determining  Mode  II  delamination 
toughness.  Results  showed  that  the  4ENF  test  yields  toughness  similar  to  that 
obtained  by  ENF  testing.  The  ENF  tests  performed  by  many  researchers  [29-32] 
showed  that  stitching  indeed  increases  the  delamination  toughness  of 
composites.  Flowever,  because  of  the  limitations  of  the  ENF  setup,  a fully 
developed  bridging  zone  is  never  achieved.  Thus  the  instantaneous  strain 
energy  release  rate  will  keep  changing  as  the  crack  propagates.  Chen,  Ifju  and 
Sankar  [33]  used  a new  setup  (cantilever  beam  setup)  to  obtain  much  longer 
crack  length.  They  incorporated  high  sensitivity  moire  interferometry  to 
determine  the  crack  tip  and  the  relative  displacement  of  the  two  crack  surfaces. 
Chapter  4 gives  the  detailed  descriptions  of  the  new  setup  and  test  results  and 
analysis. 

Ridards  and  Korjakin  [34]  used  the  traditional  mixed-mode  setup  to  test 
the  fracture  toughness  of  unstitched  laminated  composite.  The  ratio  of  Mode  I 
and  Mode  II  was  fixed  at  1.33  according  to  linear  beam  analysis.  Flowever  the 
ratio  is  not  constant  as  the  crack  propagates.  Reeder  and  Crew  [35-36] 
proposed  a new  mixed-mode  bending  method  for  delamination  testing.  In  this 
setup,  a single  applied  load  simultaneously  produces  a Mode  I and  Mode  II 
bending  load  on  the  specimen.  This  test  offers  numerous  ratios  of  Mode  I and 
Mode  II  and  has  several  advantages  over  the  traditional  mixed-mode  test. 
Shivakumar  et  al.  [37]  used  a Modified  Mixed  Mode  Flexure  Specimen  for  the 
Mixed  Mode  I/ll  tests  and  obtained  the  total  critical  strain  energy  rate  Gc.  Then 
Mode  I Gic  and  Mode  II  GnC  were  separated  by  linear  beam  theory  or  by  FE 
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analysis,  which  yields  slightly  higher  mixed  mode  ratio  and  some  dependence  on 
crack  length.  Results  showed  some  differences  between  the  values  of  the  critical 
release  rates  for  crack  initiation  and  aftermath.  Initial  Gc  was  lower  than  the 
average  G.  Chen  et  al.  [38]  modified  the  mixed-mode  bending  test  apparatus 
and  eliminated  the  error  caused  by  the  weight  of  the  lever.  Cox  and  Massabo 
[39]  investigated  the  mixed  mode  delamination  behavior  of  through-thickness 
carbon-epoxy  laminates  by  using  two  different  test  specimens  (T-stiffener  and 
mixed-mode  bending  specimen).  They  found  that  small  quantities  of  titanium  or 
carbon  z-fiber  substantially  improve  delamination  resistance  in  both  types  of 
specimen.  Reinforcement  raises  the  ultimate  strength  of  the  MMB  specimen  by 
a factor  of  three.  Although  many  researchers  [40-46]  used  different  approaches 
to  investigate  the  delamination  fracture  toughness  of  composites  and  to  obtain 
useful  information  about  the  mixed  mode  fracture  properties  of  composites,  no 
researchers  have  successfully  tested  the  mixed  mode  fracture  toughness  of  the 
composites  with  dense  translaminar  reinforcement.  Ifju,  Chen  and  Sankar  [47] 
developed  three  new  mixed  Mode  setups  that  can  be  used  to  measure  any 
conceivable  ratio  of  Mode  I and  Mode  II.  Detailed  setups  and  test  results  are 
shown  in  Chapter  5. 

Byun  et  al.  [48]  conducted  a finite  element  analysis  of  3-D  woven  double 
cantilever  beam  (DCB)  specimens  and  evaluated  the  Mode  I critical  strain  energy 
release  rate  to  investigate  the  influence  of  through-the-thickness  fibers  on  crack 
driving  force  and  crack  length.  Chen  et  al.  [49]  proposed  the  effective  critical 
strain  energy  release  rate  to  measure  Mode  I fracture  toughness  of  stitched 
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laminates  using  a finite  element  model.  Jain  and  Mai  [50]  have  analytically 
modeled  the  Mode  I and  Mode  II  delamination  toughness  of  stitched  laminated 
composites.  Sankar  and  Zhu  [51]  used  analytical  simulation  to  characterize  the 
stitch’s  effect  on  laminated  composites.  They  confirmed  that  stitching  greatly 
increases  Mode  I fracture  toughness.  Sankar  and  Dharmapuri  [52]  used  a 
simple  analytical  model  to  describe  crack  propagation  in  stitched  laminated 
Double  Cantilever  Beam  (DCB)  specimens  to  obtain  a closed  form  solution. 
Compared  to  experimental  results,  it  is  shown  that  inelastic  behavior  of  the 
stitches  plays  a significant  role  in  increasing  fracture  toughness  due  to  stitching. 
Rhee  and  Lee  [53]  used  a two-term  parameter  technique  to  calculate  the  energy 
release  rate.  Since  this  technique  uses  far-field  stress  and  displacement 
distributions  from  the  delamination  tip,  it  does  not  require  a refined  mesh  near  the 
crack  tip.  Poursartip  and  Chinatambi  [54]  presented  a method  for  experimental 
determination  of  the  local  strain  energy  rate  based  on  the  crack  opening 
displacement  (COD)  and  crack  shear  displacement  (CSD)  profiles  at  the  crack  tip 
for  the  case  of  Mode  I and  Mode  II  loading  of  composite  laminates.  It  was  found 
that  the  magnitudes  of  the  measured  COD  and  CSD  profiles  agree  well  with 
analytical  linear  elastic  fracture  mechanics  (LEFM)  predictions.  Numerous  other 
models  [55-61]  were  proposed  to  characterize  the  composites’  properties.  These 
models  make  it  convenient  to  change  one  or  more  parameters  of  composites. 
However,  their  models  only  give  one  parameter  to  characterize  the  fracture 
toughness  of  stitched  composites.  In  Chapter  6,  the  author  proposes  two 
parameters  to  describe  the  fracture  toughness  of  stitched  composites  - one  is 
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Gparent  which  is  only  related  to  the  parent  material  and  is  constant  throughout  a 
structure.  Another  is  Geff  which  is  not  only  related  to  the  parent  material  but  also 
depends  on  the  stitches  in  the  bridging  zone.  A comparison  of  these  two 
parameters  will  show  the  effectiveness  of  stitching  on  the  structure. 
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Stitching  Z-pinning 

Figure  1-1.  Two  types  of  translaminar  reinforcement:  continuous  (Stitching)  and 
discontinuous  (Z-pinning). 


13 


(A) 


■— — 

— V 



1 

(B) 


Figure  1-2.  Types  of  Stitches:  (A)  lock  stitch,  (B)  modified  lock  stitch,  and  (C) 
chain  stitch. 


CHAPTER  2 

PRELIMINARY  MODE  I TEST  AND  ANALYSIS 
Through-the-thickness  stitching  dramatically  improves  the  interlaminar 
fracture  toughness  of  laminated  composites.  Currently  the  double  cantilever 
beam  (DCB)  test  is  used  to  determine  the  Mode  I interlaminar  fracture  toughness 
of  unstitched  and  lightly  stitched  composite  laminates.  But  the  standard  DCB 
test  method  is  not  suitable  for  laminates  with  medium  to  high  density  stitching. 
Because  of  the  high  bending  moment  at  the  crack  tip  the  specimens  fail  due  to 
high  compressive  stresses  before  the  delamination  can  propagate.  To  overcome 
this  limitation,  a preliminary  test  setup,  which  induces  axial  tension  in  the 
specimen  in  addition  to  the  transverse  force  responsible  for  propagation  of 
delamination,  has  been  developed.  The  axial  tension  reduces  the  compressive 
stresses  in  the  vicinity  of  the  crack  tip  caused  by  the  large  bending  moments 
required  for  crack  propagation.  An  analysis  method  has  been  developed  to 
determine  load-deflection  and  load-energy  release  rate  relationships  in  the  new 
test.  The  new  test  method  has  been  found  to  be  successful  in  testing 
graphite/epoxy  laminates  containing  Kevlar  stitches  with  a density  of  64  stitches 
per  square  inch.  The  apparent  fracture  toughness  of  these  specimens  is  about 
45  times  that  of  unstitched  specimens. 
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2.1  Introduction 

The  standard  DCB  specimen  (Figure.  2-1)  is  fabricated  with  a Teflon  layer 
at  the  mid-plane  to  simulate  an  initial  crack.  When  loaded,  the  crack  propagates 
in  the  plane  of  delamination.  Mode  I fracture  toughness  can  be  determined  using 
one  of  the  two  following  methods.  The  load  at  which  the  delamination  starts 
propagating  can  be  used  to  compute  the  Gic  as 


G,c  ~ 


py_ 
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where  Fc  is  the  load  at  which  the  crack  propagation  occurs,  a is  the  current  crack 
length,  b is  the  width  and  El  is  the  equivalent  flexural  rigidity  of  the  specimen.  In 
the  second  method  the  specimen  is  unloaded  after  the  crack  propagates  through 
some  distance,  say  Aa.  The  area  under  the  load-deflection  diagram  represents 
the  work  done  AW  in  propagating  the  crack,  and  the  critical  energy  release  rate  is 
computed  as 

^ _AW 

(_7  is*  — 

bAa 

For  instance  Gic  of  T300/3501-6  unidirectional  unstitched  graphite/epoxy 
laminates  is  about  1.7  lb  in/in2  (298  J/m2).  It  should  be  noted  when  testing 
stitched  specimens,  unloading  can  never  be  complete  as  the  broken  stitches 
protrude  out  of  the  newly  created  delamination  surface  and  prevent  them  from 
closing.  In  that  case  the  author  assumes  that  the  unloading  is  elastic  and 
connect  the  current  point  in  the  load-deflection  curve  to  the  origin. 

Unfortunately,  the  conventional  DCB  test  is  not  suitable  for  medium  to 
high  density  stitched  composites.  The  reason  is  as  follows.  Stitches  are  very 
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strong  and  large  forces  are  required  to  break  them  or  cause  pullout.  When  the 
load  increases,  the  specimen  arms  are  subjected  to  large  bending  moments.  In 
the  end,  a specimen  arm  fails  because  of  micro-buckling  on  the  compression 
side  before  crack  propagation  as  shown  in  Figure.  2-2.  The  energy  release  rate 
in  a double  cantilever  beam  can  be  written  as  the  difference  in  strain  energy 
densities  (strain  energy  per  unit  length  of  the  beam)  behind  and  ahead  of  the 
crack  tip: 


bEI 


where  M is  the  bending  moment  just  behind  the  crack  tip.  Thus  the  bending 
moment  required  to  propagate  the  crack  varies  as  M = (G/cbEI)1/2.  Because 
stitched  laminates,  even  with  low  stitch  density,  exhibit  about  25  times  increase 
in  fracture  toughness,  the  bending  moment  required  is  5 times  that  for  unstitched 
laminates.  This  factor  is  even  higher  for  medium  to  heavily  stitched  laminates. 
Such  a high  bending  moment  induces  enormous  axial  stresses  on  the  top  and 
bottom  surfaces  of  the  sub-laminates.  Because  unidirectional  graphite/epoxy 
composite  has  very  high  tensile  strength,  the  high  tensile  stresses  induced  by  the 
bending  moment  do  not  seem  to  be  a problem.  However  failure  initiates  on  the 
compressive  side  before  the  crack  can  propagate. 

Previous  researchers,  for  example,  have  suggested  incorporation  of 
reinforcing  tabs  as  shown  in  Figure  2-3.  Aluminum  tabs  bonded  to  the  top  and 
bottom  surfaces  of  the  specimen  can  increase  the  load  carrying  capability.  In 
theory,  this  allows  the  specimen  to  resist  larger  bending  moment.  But,  in 
practice,  this  scheme  faces  another  problem  — when  loaded,  the  aluminum  tabs 
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debond  from  the  specimen  due  to  high  interlaminar  shear  stresses.  Further  they 
can  alter  the  stitch  failure  mechanism  and  the  measured  Gic  may  not  be  the  true 
fracture  toughness. 

2.2  Modified  DCB  Test 

2.2.1  New  Fixture 

To  avoid  the  aforementioned  compressive  failure,  one  must  reduce  the 
maximum  compressive  stresses.  By  adding  tensile  stresses  along  the  in-plane 
fiber  direction,  the  maximum  compressive  stresses  can  be  suppressed. 
Following  this  idea,  a preliminary  fixture  was  developed  as  shown  in  Figure  2-4. 
In  this  fixture,  heavy  aluminum  tabs  are  bonded  to  the  specimen  ends.  Kevlar 
yarns  are  wound  over  the  tabs  and  the  specimen  to  prevent  the  specimen  from 
pulling  out  of  the  tabs.  Two  rollers  attached  to  the  tabs  roll  over  a pair  of 
adjustable  rails.  As  the  specimen  is  pulled  vertically,  the  rollers  follow  the  path  of 
the  rails  opening  the  specimen  ligaments.  Thus  the  transverse  force  that  opens 
the  crack  is  proportional  to  the  tensile  force  applied  to  the  specimen.  The  tensile 
force  effectively  cancels  or  reduces  the  compressive  stresses  caused  by  the 
bending  moment  behind  the  crack  tip.  Bearings  are  used  in  the  rollers  to  reduce 
friction.  Additionally,  the  angles  of  the  rails  are  adjustable  and  hence  the  ratio 
between  the  vertical  force  and  the  horizontal  force  can  be  adjusted.  For  a given 
stitched  specimen,  the  optimal  angle  can  be  determined.  For  too  small  angle  of 
the  rail  (defined  by  a in  Figure  2-4)  there  is  too  much  transverse  force  and  the 
specimen  breaks  before  the  stitches  break.  For  too  large  a the  transverse  force 
is  small  and  the  specimen  has  to  undergo  large  axial  displacement  before  the 
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crack  can  propagate.  Thus  we  want  to  find  the  maximum  a that  will  be  able  to 
break  the  stitches  without  breaking  the  specimen. 

2.2.2  Specimen  Preparation 

The  stitched  specimens  were  made  of  24  plies  of  AS4  uniweave  graphite 
fabric  and  3501-6  epoxy  resin  using  the  RTM  process.  Stitch  yarns  (bobbin  yarn) 
were  1600  denier  Kevlar  (2790  yard/lb).  Stitch  density  of  the  specimen  was 
8x1/8  inch.  We  define  the  stitch  density  by  the  number  of  stitches  per  square 
inch  and  represent  this  density  by  the  stitching  pattern  as  (Number  of  stitches  per 
inch)x(Spacing  between  two  stitch  lines)  (e.g.,  8x1/8  inch  means  a stitch  density 
of  64  where  pitch  is  8 stitches  per  inch  and  distance  between  two  adjoining  stitch 
rows  is  1/8  inch).  Needle-stitching  yarn  used  in  all  the  cases  was  Kevlar-29.  Top 
and  bottom  plies  of  the  uniweave  preform  were  covered  by  one  layer  of  plain- 
weave  fiberglass  cloth  to  act  as  the  retainer  cloth  for  the  stitches.  A typical  test 
specimen  was  about  0.26  inch  wide,  0.16  inch  thick  and  7 inches  long.  There 
are  two  rows  of  stitches  in  each  specimen.  The  bonded  areas  were  filed  to 
increase  the  friction  and  the  Kevlar  threads  were  used  to  bond  the  aluminum  tabs 
to  the  specimen.  Screws  were  used  to  tighten  the  tab  with  the  specimen  and 
produce  very  high  pressure  between  the  tabs  and  the  surfaces  of  the  specimen 
to  prevent  slipping  of  the  tabs  from  the  bonded  area  of  the  specimen. 

2.2.3  Fixture  Verification 

Two  sets  of  unstitched  specimens  were  cut  from  the  same  unidirectional 
24  ply  graphite/epoxy  panel.  Each  set  includes  4 specimens.  One  set  of 
specimens  were  tested  by  the  standard  DCB  method  and  the  other  by  the  new 
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fixture.  Table  2.1  shows  the  comparison  between  the  two  methods.  Gic  from  the 
new  method  is  slightly  greater  than  that  from  conventional  DCB  tests.  Because 
the  new  fixture  incorporates  rails  and  bearings,  the  friction  in  the  mechanisms 
consumes  some  energy.  Currently  the  fixture  design  is  being  modified  so  that 
friction  can  be  kept  to  a minimum. 


Table  2.1 . Comparison  of  two  methods 


STANDARD  DCB  TESTS 

NEW  FIXTURE  TESTS 

Specimenl  G|C=1. 63  lb  in/in2 
Specimen2  Gic=1.81  lb  in/in2 
Specimen3  G|C=1. 57  lb  in/in2 
Specimen4  Gic=1. 67  lb  in/in2 

Average:  G|C=1.67  lb  in/in2 

Specimenl  Gic=1. 78  lb  in/in2 
Specimen2  G|C=1. 97  lb  in/in2 
Specimen3  Gic=1. 99  lb  in/in2 
Specimen4  Gic=1. 77  lb  in/in2 

Average:  Gic  =1.88  lb  in/in2 

2.2.4  Testing  of  Stitched  Specimens 

Tests  were  conducted  in  a screw-driven  universal  testing  machine.  The 
crack  propagation  was  observed  by  a CCD  camera.  A computer  was  used  to 
record  the  images  corresponding  to  various  displacements.  The  load,  crack 
extension  and  crack  opening  were  determined  corresponding  to  each 
displacement  step.  Figure  2-4  shows  when  the  grip  of  the  machine  moves  up, 
the  bearings  on  the  specimen  follow  the  rails  of  the  fixture  and  open  up  the  crack 
surfaces.  The  optimum  angle  of  the  rails  depends  on  the  mechanism  of 
particular  stitched  composites.  The  stitches  typically  break  or  pull  through  the 
composite.  The  horizontal  components  of  the  reactions  on  the  bearings  tend  to 
open  the  crack  while  the  vertical  component  provides  the  tensile  forces  to  reduce 
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the  compressive  stresses  induced  by  the  bending  moment.  Using  this  fixture,  we 
were  able  to  test  specimens  with  medium  as  well  as  very  high  stitch  densities. 

Typical  measurements  obtained  using  the  new  fixture  are  shown  in  Figure 
2-6.  During  the  test  the  load,  displacement  in  the  direction  of  applied  load,  crack 
surface  opening  displacement  and  the  crack  length  are  all  measured.  It  should 
be  mentioned  that  the  crack-surface  opening  displacement  is  measured  at  a 
reference  point  near  the  end  of  the  tabs.  The  crack  length  is  defined  by  the 
distance  between  this  reference  point  and  the  crack-tip.  As  a sample,  the  load, 
crack  length  and  crack  opening  displacement  are  plotted  against  the  axial 
displacement  for  one  of  the  tests  in  Figure  2-6.  From  the  load-displacement 
diagram  in  Figure  2-6,  one  can  note  that  initially  the  load  increases  almost 
linearly  with  the  displacement.  When  the  stitches  break,  the  crack  begins  to 
propagate  and  there  is  a sudden  drop  in  the  load.  As  the  loading  is  continued, 
the  load  goes  up  and  down  as  the  stitches  break  and  as  the  crack  propagates 
intermittently.  This  saw  tooth  behavior  of  the  load-deflection  curve  is  due  to  the 
stitches  breaking  and  new  stitches  offering  resistance  to  crack  propagation.  The 
area  under  the  load  deflection  curve  represents  the  energy  used  in  breaking  the 
stitches  and  propagating  the  delamination. 
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2.3  Analysis 

2.3.1  Governing  Equations 

The  Gic  can  be  determined  using  one  of  two  methods.  In  the  first  method, 
the  area  under  the  load-deflection  diagram  is  divided  by  the  crack  extension 
area: 


G.C  ~ 


AW 
b ■ A a 


(1) 


where  AW  is  the  work  of  fracture,  A a is  the  crack  extension  and  b is  the  width  of 
the  specimen.  This  method  relies  completely  on  experiments.  As  mentioned 
earlier  the  energy  released  during  unloading  is  obtained  by  joining  the  current 
point  in  the  load  deflection  curve  to  the  origin. 

The  second  method  involves  an  elastic  beam  analysis  of  the  specimen.  A 
large  deflection  nonlinear  analysis  is  performed  to  derive  the  load-deflection 
relationship  of  the  test  specimen.  Then  from  the  analysis  results,  an  expression 
for  the  strain  energy  in  the  beam  can  be  derived.  The  critical  energy  release  rate 
at  any  given  load  (or  deflection)  then  can  be  obtained  as: 

dU 


G]C  — 


b ■ da 


(2) 


where  U is  the  strain  energy  in  the  beam  as  a function  of  crack  length  a, 
load/deflection  and  other  beam  properties.  The  analysis  procedure  is  described 
below. 


Due  to  symmetry  one  of  the  ligaments  of  the  DCB  specimen  (Figure  2-4) 
are  modeled.  Timoshenko  beam  theory,  which  accounts  for  transverse  shear 
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deformation,  is  used.  The  two  governing  equations  in  the  variables  w and  y, 
transverse  deflection  and  rotation  respectively,  are: 


A55  ('f'  + ~)  = Q 
OX 

D^-  = M + P(w0-w)-Q-(a-x) 
dx 


(3) 


where  A55  is  the  transverse  shear  stiffness,  D is  the  flexural  stiffness,  P is  the 
axial  load,  and  Q and  M are  the  shear  force  and  bending  moment  at  a given 
cross  section,  and  w0  is  the  transverse  deflection  at  the  aforementioned 
reference  point.  It  may  be  noted  that  the  term  P(w0-w)  accounts  for  the  stress 
stiffening  of  the  beam  due  to  the  axial  force  P.  The  solution  to  these  equations 
can  be  derived  as  equation  (4). 


w = C,  • sinh(T  x)  + C2  ■ cosh(A  x)  + 
Q x | (M  + P ■ w0  - Q ■ a) 

P P 


Q(- 


ct  = 


*55 


4) 


C M + P-w0-Q-a ) 
P 


A = 


(4) 


Once  expressions  for  w(x)  and  y(x)  are  obtained  the  force  and  moment 
resultants  can  be  easily  derived  using  the  beam  constitutive  relations.  From  the 
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results  the  strain  energy  in  the  beam  can  be  calculated  using  numerical 
integration  of  the  expression  (5). 

Once  expressions  for  w(x)  and  v|/(x)  are  obtained  the  force  and  moment 
resultants  can  be  easily  derived  using  the  beam  constitutive  relations. 

From  the  results  the  strain  energy  in  the  beam  can  be  calculated  using 
numerical  integration  of  the  following  expression: 


P 2 

2 bhE 


Q 


2 A 


2bA 


dx 


55  y 


(5) 


Then  the  critical  energy  release  rate  in  Equation  (2)  can  be  obtained  by 
numerical  differentiation  of  U(a)  with  respect  to  the  crack  length  a 

2.3.2  Verification  of  the  Model 

Figure  2-7  illustrates  the  comparison  of  the  analytical  solution  with  the  test 
data.  The  dashed  line  represents  the  crack  extension  vs.  crack  opening  data 
from  a test.  We  used  the  loads  and  crack  extension  from  the  test  data  to 
calculate  the  crack  surface  opening  which  is  plotted  as  a solid  line  in  Figure  2-7. 
The  relative  error  between  the  measured  opening  and  the  calculated  opening  is 
below  10  percent.  Because  our  analysis  results  are  close  to  the  experimental 
data,  we  use  this  method  with  confidence  to  obtain  the  instantaneous  critical 
strain  energy  release  rate  Gic- 

2.3.3  Instantaneous  Gin 

In  this  method  the  analysis  described  above  is  used  to  obtain  two  sets  of 
load-deflection  curves,  one  for  a crack  length  a and  another  for  a+Aa,  where  Aa 
is  a small  but  arbitrary  increment.  A set  of  sample  load-deflection  curves  is 
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shown  in  Figure  2-8.  Then  from  the  two  curves  the  critical  energy  release  rate 
can  be  obtained  numerically  using  Equation  (2).  The  variation  of  G|  as  a function 
of  the  load  for  a given  crack  length  is  shown  in  Figure  2-9.  The  value  of  G| 
corresponding  to  the  load  at  which  the  crack  began  to  propagate  is  the  fracture 
toughness  of  the  stitched  specimen  at  that  instant.  This  procedure  can  be 
repeated  for  various  crack  lengths  and  the  fracture  toughness  can  be  determined 
as  a function  of  the  crack  length.  For  this  stitched  specimen,  Figure  2-10  shows 
the  fracture  toughness  remains  almost  constant  with  respect  to  the  crack  length 
and  is  approximately  equal  to  80in  lb/in2  (14,000  J/m2).  By  comparing  to  the  G|C 
of  unstitched  specimens  (Table  2-1),  one  can  note  that  the  Mode  I fracture 
toughness  increases  by  more  than  45  times  due  to  stitching. 

2.4  Conclusions 

A new  preliminary  DCB  test  has  been  developed  for  Mode  I fracture 
testing  of  laminated  composites  containing  a high  density  of  through-the- 
thickness  stitches.  Conventional  DCB  tests  are  not  suitable  for  stitched 
specimens,  since  the  specimen  breaks  just  behind  the  crack-tip  due  to  high 
bending  stresses,  especially  compressive  stresses.  In  the  present  method  a 
tensile  stress  proportional  to  the  transverse  force  is  induced  which  reduces  the 
amount  of  compressive  stresses.  An  analysis  method  is  developed  to  compute 
the  Mode  I fracture  toughness  from  the  load-deflection  relation.  The  method  has 
been  validated  by  testing  unstitched  specimens  for  which  the  fracture  toughness 
is  known,  and  then  used  to  test  stitched  graphite/epoxy  composites.  The  Gic 
values  for  unstitched  specimens  were  slightly  greater  in  the  present  method.  This 
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is  due  to  frictional  dissipation  in  the  fixture.  The  apparent  fracture  toughness  of 
specimens  tested  was  found  to  be  about  45  times  that  of  unstitched  specimens. 
Efforts  are  underway  to  modify  the  fixture  to  reduce  friction  and  improve 
repeatability. 
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Figure  2-1.  Standard  DCB  test  method 
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Figure  2-2.  Standard  DCB  test  for  high  density  stitched  composites 
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Figure  2-4.  Modified  DCB  test 
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Figure  2-5.  Progressive  loading  of  the  preliminary  fixture 
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Figure  2-6.  Axial  displacement  vs.  load,  crack  extension  and  opening 
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Figure  2-7.  Crack  extension  vs.  opening 
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Figure.  2-8.  Load  vs.  axial  displacement  under  two  crack  lengths  - a and  a’ 
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Figure  2-9.  Critical  strain  energy  release  rate  vs.  load 
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Figure. 2-10.  Critical  strain  energy  release  rate  (G|C)  vs.  crack  length 


CHAPTER  3 

FINAL  DESIGN  OF  MODE  I TEST  FIXTURE 


Previous  chapter  introduced  the  preliminary  Mode  I test  for  stitched 
composites.  It  has  been  used  to  test  the  specimens  with  stitching  density  up  to 
64  stitches  per  square  inch.  However,  this  fixture  has  some  disadvantages,  for 
instance,  the  poor  gripping  method  and  the  friction  between  the  railings  and 
bearings.  After  trial  and  error,  the  author  developed  a new  fixture  which 
overcomes  the  deficiencies  of  the  prototype  design.  It  can  be  very  easily  gripped 
to  a specimen  and  can  be  adjustable  to  different  thickness  of  specimens.  Most 
importantly,  there  is  no  friction  involved  by  using  this  new  fixture.  This  fixture  has 
the  potential  to  be  the  standard  test  method  for  laminated  composites  with 
translaminar  reinforcement.  Experiments  were  conducted  using  carbon/epoxy 
specimens  containing  6.2  stitches  per  square  centimeter  (40  stitches  per  square 
inch).  The  nonlinear  differential  equations  of  equilibrium  of  the  new  specimen 
are  solved  using  an  iterative  procedure  to  obtain  the  strain  energy  release  rate  as 
a function  of  load  and  crack  length.  Results  include  Mode  I fracture  toughness, 
crack  tip  bending  moment,  transverse  deflection  and  slope  as  a function  of  crack 
length.  It  is  found  that  the  apparent  fracture  toughness  of  the  specimens  tested 
remains  constant  as  the  stitches  break  and  the  crack  propagates,  and  is  about 
sixty  times  that  of  unstitched  specimens. 
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3.1  Description  of  the  New  Test  Fixture 
The  standard  DCB  specimen  is  depicted  in  Figure  2-1.  The  energy 

M2 

release  rate  in  the  DCB  specimen  is  given  by , where  M = Qa  is  the  bending 

bEI 


moment  in  one  of  the  sub  laminates  at  the  crack  tip,  Q is  the  transverse  force,  a 
is  the  crack  length,  b is  the  width  of  the  specimen  and  El  is  the  flexural  rigidity  of 
the  sublaminates.  As  the  fracture  toughness  of  the  material  increases,  larger 
bending  moments  are  required  to  propagate  the  crack.  The  bending  moments 
also  produce  large  compressive  stresses  on  the  outer  surfaces  of  the  specimen 
as  indicated  in  Figure  2-2.  Although  fiber  composites  are  strong  in  tension,  their 
compressive  strength  is  limited  by  fiber  micro-buckling,  which  causes  the 
specimens  to  fail  before  the  crack  could  propagate.  One  of  the  ways  of  reducing 
the  compressive  stresses  is  to  apply  a uniform  tension  in  both  arms  of  the  DCB 
specimen.  It  should  be  noted  that  this  tensile  stress  being  equal  in  both  arms 
does  not  contribute  to  the  energy  release  rate,  and  Mode  I crack  propagation  is 
maintained.  Following  this  idea  a new  fixture  has  been  developed  as  illustrated 
in  Figure  3-1.  This  fixture  can  supply  forces  along  two  directions.  The 
transverse  component  of  the  force  is  for  crack  opening  and  the  axial  component 
is  for  applying  the  tensile  force.  The  fixture  has  been  designed  such  that  the 
initial  angles  of  the  arms  are  adjustable  and  for  a given  stitched  specimen,  the 
optimal  angle  can  be  determined. 

A challenge  in  the  new  fixture  is  the  design  of  the  specimen  grips  that  can 
transmit  both  axial  and  transverse  forces.  The  conventional  method  of  bonding 
the  tabs  to  the  specimen  does  not  work  as  the  bonded  tabs  cannot  withstand  the 
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large  tensile  and  shear  stresses  caused  by  the  larger  load  required  to  propagate 
the  crack.  A notch  in  the  form  of  a circular  arc  was  machined  in  the  specimen 
ends  as  shown  in  Figure  3-2.  The  notch  surface  forms  a 45-degree  angle  to  the 
specimen  surface.  A pair  of  adjustable  grips  that  match  the  notch  profile  in  the 
specimen  were  machined  out  of  steel  (Figure  3-3.).  The  final  assembly  of  the 
specimen  and  the  grips  is  shown  in  Figure  3-4.  Figure  3-5  shows  a series  of 
pictures  taken  from  one  of  the  tests.  The  fixture  unzips  the  stitches  one  by  one 
and  the  crack  propagation  is  found  to  be  stable. 


A schematic  of  the  new  DCB  test  set-up  and  the  free-body  diagram  of  AB 
are  shown  in  Fig.  3.6.  As  a row  of  stitches  breaks,  the  crack  front  reaches  the 
location  of  the  next  row  of  the  stitches.  The  current  crack  tip  is  located  at  Point  B, 
and  the  portion  of  the  specimen  in  the  grip  (AC  in  Fig.  3.6)  is  considered  rigid.  All 
of  the  forces  and  the  bending  moment  at  A can  be  expressed  as: 


= FL3  sin(<9,  -02)~  FL4  cos(0]  -02) 

The  angles  Oi  and  02  are  indicated  in  Figure  3-6.  The  angle  02  is  the  slope 


3.2  Analysis  of  the  Specimen 


fa=f 

Ma  = F(L3  cos  9 2 - L4  sin  02 ) sin  0X  - F(L3  sin  02  + L4  cos  02 ) cos  0X 


(1) 


at  point  A,  i.e.,  02  - -----  X=1  . The  angle  0i  can  be  calculated  using  the 

/7v  5 


dx 


following  geometrical  relations: 


=Ao  +L3  *sin02  +vA,  vA=v(x)x=Ls, 


(2) 


where  the  distance  A0  is  indicated  in  Figure  3-1. 
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One  only  needs  to  analyze  the  portion  AB  as  shown  in  Figure  3-6.  The 

expression  for  the  bending  moment  M(x)  can  be  derived  as 

M(x ) = M A + F(L5  - x ) sin  9X  - F • [v^  - v(x)]  • cos  Gx 

= FL3  sin(0j  -02)~  FL 4 cos(0,  -02)+  (3) 

F(Z5  - x) sin 6X  -F-[va  - v(x)]  • cos 

where  va  and  v(x)  are  transverse  deflections  at  cross  section  A and  an  arbitrary 
cross  section  at  x.  It  must  be  mentioned  that  the  calculation  of  bending  moment 
M(x)  in  Equation  (3)  is  exact,  and  the  contribution  of  the  axial  force  to  the 
bending  moment  due  to  large  deflection  is  included.  The  governing  differential 
equation  for  the  beam  is: 

1 d2v(x)  _ 1 d0  _ M(x)  . 

(\  + (dv/dx)2f/2  dx 2 ~ (\  + 02y,2~fa~  El 

where  El  is  the  equivalent  flexural  rigidity  of  one  arm  of  the  specimen  and  the 

expression  on  the  left  hand  side  of  the  equation  is  exact  for  the  beam  curvature. 

At  this  point  we  make  the  usual  approximation  for  the  curvature: 

d2v(x ) _d0  _ M(x) 
dx2  El 

From  the  experiments  (Figure  3-5)  we  noted  that  the  slope  of  the  beam  is 
higher  towards  the  tip  of  the  DCB  specimen  where  the  bending  moments  are 
small.  On  the  other  hand,  towards  the  crack  tip,  where  the  bending  moment  is 
large,  the  slopes  are  indeed  small.  Hence  we  are  justified  in  neglecting  the 
higher  order  term  in  the  denominator  for  the  curvature  and  this  further  simplifies 


the  calculations. 
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The  equivalent  flexural  rigidity  El  was  measured  from  a three-point 
bending  test  in  which  the  maximum  deflection  8 at  the  center  was  related  to  the 
applied  force  P by  8 = PL3/48EI  where  L is  the  span  of  the  loading  fixture.  The 
bending  moment  M(x)  is  not  only  a function  of  x,  but  also  a function  of  0,  and  02 . 
An  iterative  numerical  method  was  used  to  solve  Equation  (3).  The  load  F and 
the  crack  length  L5  are  obtained  from  the  test  data.  First  the  initial  values  of  91t 
02  and  vA  are  calculated  for  the  initial  position  of  the  specimen  and  Equation  (5) 
was  solved  for  v(x).  After  the  first  iteration,  new  values  of  Qi,  92  and  vA  were 
calculated  and  used  for  the  second  round  of  the  iteration.  As  the  relative  errors 
between  successive  sets  of  values  of  91t  92  and  vA  become  less  than  an 
allowable  minimum,  say  2%,  the  iterations  are  stopped,  and  the  final  numerical 
values  are  taken  as  the  solution.  Following  are  steps  used  in  the  iterative 
process: 

Stepl : import  F (load),  L5  (crack  length)  from  the  test. 

Step2:  use  the  initial  values  of  9-,,  92  and  vA  to  solve  v(x). 

Step3:  from  v(x),  obtain  new  91t  92  and  vA. 

Step4:  calculate  the  relative  errors  between  new  9i,  92  and  vA  and  previous  9i,  92 
and  vA.  If  the  errors  are  allowable,  stop  and  output  new  91t  92  and  vA  as  the  final 
values  of  9i,  92  and  vA.  Otherwise,  go  to  step2  by  using  new  9h  92  and  vA  as 
input  to  calculate  v(x). 
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From  the  results,  the  bending  moment  and  shear  force  at  the  section  just 
behind  the  crack  tip  can  be  calculated  for  a given  load  and  crack  length.  The 
energy  release  rate  is  obtained  from  the  formula: 


GIC  =- 
,c  b 


kEI  j 


(6) 


where  b is  the  width  of  the  specimen  and  M is  the  crack  tip  bending  moment  at 
the  instant  of  crack  propagation.  It  should  be  noted  that  in  order  to  use  the 
aforementioned  iterative  method,  we  need  to  measure  the  crack  length 
corresponding  to  the  peak  load  in  the  tests. 

It  must  be  mentioned  that  the  aforementioned  approach  ignores  the  crack 
bridging  effects  caused  by  the  unbroken  stitches  in  the  wake  of  the  crack  tip. 
Such  crack  bridging  effects  in  stitched  composites  have  been  considered  in  great 
detail  by  Massabo  and  Cox  [26],  Sankar  and  Dharmapuri  [52]  also  analytically 
studied  the  effects  of  bridging  in  stitched  composite  beams.  In  the  present  study 
there  are  two  reasons  for  ignoring  the  bridging  effects  (small  scale  bridging). 
First,  the  size  of  the  bridging  zone  is  small  compared  to  the  length  of  the  crack. 
Typically  it  was  about  one  stitch  spacing,  that  is,  at  any  given  time  not  more  than 
two  unbroken  stitches  were  found  behind  the  crack  tip.  Second,  this  bridging 
zone  was  found  to  be  a constant  as  the  crack  tip  moved  forward.  That  is,  there 
was  no  crack  resistance  effect  and  the  fracture  toughness  remains  constant  with 
the  crack  length.  Flence  the  present  approach  can  be  considered  as  a method  of 
determining  the  apparent  fracture  toughness  of  the  stitched  composites  where 
the  micromechanics  of  stitch  failure  are  not  considered  in  detail. 
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3.3  Results  and  Discussion 

There  are  two  methods  of  determining  G/c  from  the  DCB  tests.  In  the  area 
method  the  work  done  by  the  external  forces  on  the  specimen  is  calculated  from 
the  area  enclosed  by  the  load-deflection  diagram.  This  work  is  divided  by  the 
area  of  crack  extension  to  obtain  the  average  G/c.  We  did  not  use  this  method 
as  the  specimen  could  not  be  unloaded  completely  because  of  the  protruding 
broken  stitches  on  the  delamination  surfaces.  These  stitches  prevented 
complete  closure  of  the  two  ligaments  of  the  DCB  specimen.  Hence  it  was 


decided  to  use  the  formula  G,c 


M2 

bEI 


where  M is  the  crack  tip  bending  moment 


at  the  instant  of  crack  propagation.  The  value  of  G/c  calculated  using  this 
method  can  also  be  thought  of  as  the  instantaneous  fracture  toughness  at  the 
time  of  crack  propagation.  However  we  need  to  know  the  bending  moment  at  the 
crack  tip  from  the  load,  crack  length  and  other  geometrical  information.  This  is 
much  simpler  in  conventional  DCB  specimens  and  can  be  obtained  as  the 
product  of  the  load  and  the  crack  length.  In  the  current  specimen  we  need  to 
perform  an  analysis  to  obtain  the  bending  moment.  The  main  difference  is  the 
effect  of  the  axial  force  on  the  bending  moment,  which  now  is  a function  of  the 
transverse  deflection  of  the  specimen.  This  necessitates  an  iterative  procedure 
to  solve  for  the  deflection  of  the  beam.  The  analysis  procedure  is  described  in  the 
following. 

The  specimens  were  cut  from  a panel  which  is  made  of  4 stacks  of 
Saertex@  textile.  Each  stack  is  made  of  7 layers  with  fiber  orientations  given  by 
[45°/-45o/O0/9O0/Oo/-45°/450].  The  linear  density  of  the  Kevlar  stitch  is  1600 
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denier.  The  specimen  is  19.1  cm  (7.5  inches)  long  and  2.0  cm  (0.8  inch)  wide.  A 
sample  load-deflection  diagram  is  shown  in  Figure  3-7.  As  the  stitches  start 
breaking,  one  can  note  the  fluctuations  in  the  load  and  the  saw-tooth  like  pattern 
of  the  load-deflection  diagram.  Each  peak  in  the  load  corresponds  to  the 
breaking  of  one  row  of  stitches.  In  the  beginning  the  load  to  break  the  stitches 
increases  (stitches  1 through  10),  then  a steady  state  is  reached  and  the  stitches 
break  at  a constant  peak  load.  The  twenty  peak  loads  shown  in  Figure  3-7 
match  the  twenty  rows  of  broken  stitches.  After  carefully  observing  the 
specimen,  we  found  that  as  one  row  of  stitches  broke,  the  crack  front  reached 
the  next  row  of  stitches.  After  measuring  the  distance  between  the  initial  position 
of  the  crack  tip  and  the  location  of  (n+1)th  row  of  stitches,  we  can  obtain  the  crack 
extension  corresponding  to  nth  peak  load.  This  peak  load  and  the  corresponding 
crack  length  (L5)  were  used  as  the  input  data  in  Equation  (3) 

The  dimensions  of  the  fixture,  the  specimen  and  the  gripping  area  for  a 
particular  specimen  (Figure  3-1  and  Figure  3-6)  were  as  follows:  L?= 190  mm, 
L2= 370  mm,  L3=  57.8  mm,  L#=11  mm,  Z_5=  15.6  mm  (the  crack  length 
corresponding  to  the  first  peak  load  of  Pi= 3176  N),  b=1 9.05mm  (width  of  the 
specimen),  /?=2.98mm  (the  half  of  the  thickness  of  the  specimen).  The  initial 
values  of  9i,  02and  vA  (the  specimen  in  the  initial  position)  were  28.58°,  0°  and  0, 
respectively.  The  equivalent  Young’s  modulus  of  the  specimen  is  equal  to 
44.42GPa.  After  seven  iterations,  0?,  92  and  1/4  were  equal  to  26.9°,  8.8°  and 
1.36  mm,  respectively.  The  relative  errors  of  91:  92  and  vA  compared  to  the 
values  of  sixth  iteration  were  1.8%,  1.7%  and  1.2%,  respectively. 
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By  using  the  experimentally  measured  critical  load  and  crack  length  AB, 
the  bending  moment,  transverse  displacement  and  slope  along  the  AB  can  be 
determined  by  using  the  analytical  method.  For  instance,  we  selected  the  tenth 
peak  load  (load=5420  N)  and  measured  the  location  of  the  eleventh  row  of 
stitches  to  obtain  the  crack  length  ( >48=44.2  mm).  Figures  3-8  through  3-10 
show  the  bending  moment,  the  deflection  and  the  slope  along  the  specimen 
length  AB.  One  can  note  that  the  crack  tip  suffers  the  largest  bending  moment. 
If  no  axial  forces  were  applied,  the  specimen  would  undergo  micro-buckling  in 
the  vicinity  of  the  crack-tip  for  this  amount  of  the  bending  moment.  Flowever  an 
axial  force  of  2710  N significantly  reduces  the  compressive  stresses  and 
prevents  micro  buckling  in  the  vicinity  of  the  crack  tip.  The  values  of  the  bending 
moment  can  be  used  to  calculate  the  critical  strain  energy  release  rate  Gic  as 
shown  in  Equation  (6). 

In  Figures  3-11  through  3-13  the  crack-tip  bending  moment,  shear  force 
and  Gic  are  presented  as  a function  of  crack  length  for  three  different  specimens. 
The  average  values  for  each  specimen  are  shown  in  respective  figures.  The 
value  of  G|C  for  unstitched  specimen  was  found  to  be  about  300  J/m2  in  a 
previous  study  [1],  Thus  one  can  note  that  stitching  has  increased  the  apparent 
fracture  toughness  by  a factor  of  about  60. 

3.4  Conclusion 

A new  test  method  has  been  developed  for  measuring  the  Mode  I fracture 
toughness  of  laminated  composites  containing  medium  to  high  density  stitches. 
The  fixture  is  designed  such  that  a proportional  axial  force  is  applied  to  the 
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specimen  as  the  transverse  forces  on  the  DCB  specimen  increases.  This 
prevents  compressive  failure  of  the  specimen,  and  the  crack  propagates  by 
breaking  the  stitches.  The  nonlinear  differential  equations  of  equilibrium  are 
solved  using  an  iterative  process  in  order  to  determine  the  bending  moment  and 
shear  force  resultants  for  a given  load  and  crack  length.  From  the  force 
resultants  the  critical  energy  release  rate  can  be  computed.  The  fracture 
toughness  of  the  stitched  graphite/epoxy  specimens  is  found  to  be  about  60 
times  that  of  unstitched  specimens. 
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Figure  3-1 . New  DCB  test  setup 


Figure  3-2.  Notches  in  the  specimen 
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Figure  3-3.  Schematic  of  the  specimen  grips.  Note  that  the  specimen  grips 
match  the  notches  in  the  specimen  shown  in  Figure  3-2 
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Figure  3-4.  View  of  the  assembled  specimen 
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Figure  3-5.  Series  of  pictures  taken  during  the  new  DCB  test 
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Figure  3-6.  Schematic  of  the  new  DCB  test  setup  and  free  body  diagram  of 
portion  AB  of  a specimen 
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Figure  3-7.  Load  vs.  displacement  for  a DCB  specimen  — the  twenty  peak  loads 
match  twenty  rows  of  broken  stitches 
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Figure  3-8.  Bending  moment  distribution  along  one  arm  of  the  DCB  specimen  for 
a crack  length  of  44.2  mm  and  load  of  5420  N 


Transverse  deflection  (mm) 
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Distance  from  the  crack  tip  (mm) 


Figure  3-9.  Deflection  curve  of  one  arm  of  the  DCB  specimen  at  a load  of  5420 
N and  crack  length  of  44.2  mm 
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Figure  3-10.  Slope  vs.  distance  from  the  crack  tip  for  crack  length=44.2  mm  and 
load=5420  N 
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Figure  3-11.  Crack-tip  bending  moment  variation  as  the  crack  propagates. 
Results  for  three  specimens  are  shown 
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Figure  3-12.  Crack-tip  shear  forces  vs.  crack  length  for  three  specimens 
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Figure  3-13.  Instantaneous  Gic  of  three  specimens  as  a function  of  crack  length 


CHAPTER  4 

MODE  II  FRACTURE  TOUGHNESS  OF  STITCHED  COMPOSITES 

In  this  chapter,  two  experimental  approaches  were  used  to  quantify  the 
Mode  II  fracture  toughness  of  stitched  composites.  Currently  end  notched 
flexure  (ENF)  specimens  are  used  to  determine  the  effect  of  stitching  on  the 
Mode  II  fracture  toughness  of  laminated  composites.  In  ENF  tests  the 
delamination  does  not  open,  but  slides.  This  makes  it  difficult  to  locate  the  crack- 
tip  as  the  crack  propagates  and  hence  to  estimate  the  length  of  crack 
propagation.  This  uncertainty  introduces  significant  errors  in  the  estimated  value 
of  the  Mode  II  fracture  toughness  of  stitched  composite  laminates.  In  the  present 
study  high  sensitive  moire  interferometry  has  been  used  to  obtain  the  full-field 
displacements  as  well  as  the  crack  tip  location  in  ENF  specimens.  By  using 
linear  elastic  fracture  mechanics  for  orthotropic  material,  the  stress  intensity 
factor  can  be  determined  from  the  relative  crack  surface  displacements  behind 
the  crack  tip.  Due  to  high  sensitivity  and  high  resolution  of  this  technique,  the 
relative  displacement  and  the  crack  tip  can  be  determined  precisely. 

4.1  Introduction 

In  order  to  fully  understand  the  complex  behavior  of  laminated  composites 
with  translaminar  reinforcements,  full-field  experimental  displacement  analysis 
was  performed  on  the  microstructural  scale  of  the  materials  tested  in  this  study. 
To  accomplish  this,  moire  interferometry  (Post,  Hahn  and  Ifju  1994)  [24]  was 
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used  to  document  the  mechanism  of  translaminar  reinforcement,  and  their  effects 
on  damage  tolerance.  Moire  interferometry  is  a laser  based  optical  experimental 
stress  analysis  method  that  is  capable  of  revealing  displacements  on  the  sub 
wavelength  scale.  The  technique  of  moire  interferometry  is  characterized  by  high 
spatial  resolution,  full-field  information  and  ability  to  determine  both  normal  and 
shear  strain.  By  measuring  the  relative  difference  of  the  fringe  orders  between 
two  points,  the  relative  movement  or  slip  between  crack  surfaces  can  be 
determined  precisely.  By  plotting  the  displacement  difference  along  the  entire 
wake,  the  effect  of  stitching  could  be  understood.  Since  moire  has  very  high 
spatial  resolution,  the  location  of  the  crack  tip  can  be  determined  very  accurately. 

4.2  End-Notched  Flexure  Test 

Two  setups,  namely,  three-point  bending  setup  and  cantilever  beam  setup 
are  used  to  determine  the  Mode  II  fracture  toughness  of  stitched  composites. 
For  the  three-point  bending  setup,  in  displacement  control  for  unstitched 
composites,  when  0 < a < 0.347L,  the  crack  is  unstable,  but  when  0.347  < a < 
0.5L,  the  crack  is  stable  where  a is  crack  length  and  L is  total  span  of  the  three- 
point  bending  fixture.  For  three-point  bending  and  cantilever  beam  setups,  the 
crack  is  unstable  for  unstitched  composites  in  load  control.  Hence,  load  control 
cannot  be  used  to  obtain  stable  crack  propagation  of  unstitched  composites  in 
both  setups.  However,  stability  for  unstitched  specimens  doesn’t  seem 
applicable  for  stitched  specimen.  Our  tests  show  both  setups  give  quite  stable 
crack  propagation  for  stitched  composites. 
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4.2.1  Three-Point  Bending  Setup 

Figure  4-1  shows  the  standard  ENF  Mode  II  test  setup.  The  specimens 
were  made  of  AS4-3501  and  fabricated  with  a Teflon  layer  at  the  mid-plane  to 
simulate  an  initial  crack.  A grating,  required  for  moire  interferometry,  was 
replicated  on  the  edge  of  the  composite  specimen.  The  specimens  were  cut 
from  a panel  which  was  made  from  4 stacks  of  graphite/epoxy  prepreg.  Each 
stack  has  7 plies  which  are  oriented  at  [450/-450/00/900/00/-45°/450].  The  stitch 
density  of  the  specimen  was  6.2  stitches  per  square  centimeter  (40  stitches  per 
square  inch)  that  is  the  spacing  between  two  rows  of  stitches  was  0.51 
centimeter  (0.2  inch)  and  there  were  3.15  stitches  per  centimeter  (8  stitches  per 
inch)  in  a row.  The  specimen  was  12.7  centimeter  (5  inches)  long,  2.54 
centimeter  (1  inch)  wide  and  0.59  centimeter  (0.230  inch)  thick.  The  span 
between  two  supporting  points  was  10.2  centimeter  (4  inches)  and  the  initial 
crack  length,  or  the  distance  from  the  supporting  point  to  the  crack  tip,  was  2.54 
centimeter  (1  inch). 

4.2.2  Cantilever  Beam  Setup 

The  disadvantage  of  the  three-point  bending  setup  is  the  limitation  of  the 
length  available  for  crack  propagation.  The  maximum  crack  length  of  the  three 
point  bending  setup  was  around  2.5  cm  (1  inch)  therefore  it  was  not  possible  to 
obtain  a fully  developed  bridging  zone.  In  order  to  fully  understand  the  effect  of 
stitches,  the  crack  should  be  allowed  to  propagate  as  long  as  possible.  The 
farther  the  crack  propagates,  the  more  the  stitches  become  involved.  In  order  to 
accomplish  this,  a new  test  setup  has  been  developed  as  shown  in  Figure  4-2. 
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Using  this  setup,  one  can  obtain  more  than  5 centimeters  (2  inches)  of  the  crack 
length  which  is  equivalent  up  to  sixteen  rows  of  stitches.  The  length  of  a 
specimen  for  this  setup  was  5 centimeter  (2  inch)  longer  than  that  of  the  Three- 
Point  setup.  Because  of  large  deflection  and  rotation  of  the  specimen,  it’s  very 
difficult  to  apply  high  sensitivity  moire  interferometry  technique  for  this  setup.  As 
the  sensitivity  of  this  method  is  so  high,  small  amounts  of  deflection  and  rotation 
will  cause  the  fringe  pattern  to  become  too  dense  to  photograph.  After  trial  and 
error,  the  author  found  that  in  order  to  successfully  run  this  test,  a PEM-I 
(portable  experimental  moire  interferometer)  should  have  at  least  three 
adjustments:  (1)  it  should  be  able  to  rotate  up  to  13  degrees  to  compensate  for 
the  rotation  of  a specimen,  (2)  it  should  be  able  to  move  up  nearly  2.5  centimeter 
(1  inch)  to  make  sure  the  two  beams  cover  the  intended  area  after  a specimen  is 
deformed,  (3)  it  should  be  able  to  move  forward  and  backward  at  least  2.5 
centimeter  (1  inch)  to  trace  the  crack  tip  after  crack  propagation. 

Figure  4-3  shows  moire  interferometry  fringe  patterns  on  the  edge  of  an 
ENF  specimen.  Both  horizontal  (U)  and  vertical  (V)  displacement  fields  are 
presented.  By  measuring  the  fringe  order  difference,  the  relative  displacement  of 
the  two  crack  surfaces  could  be  determined  as  a function  of  position  along  the 
crack. 

4.3  Methods  of  Analysis  and  Results 

Two  different  methods  were  used  to  estimate  the  Mode  II  fracture 
toughness  from  the  experimental  results.  The  first  one  is  the  standard  area 
method  used  in  conjunction  with  ENF  tests,  and  it  yields  the  average  Mode  II 
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critical  energy  release  rate.  The  second  method  is  based  on  the  crack  surface 
displacements  measured  using  high  sensitivity  moire  interferometry.  This 
method  yields  the  Mode  II  critical  stress  intensity  factor.  While  the  area  method 
gives  the  average  Gnc  over  the  length  of  crack  propagation,  the  moire  method 
results  in  instantaneous  Knc  for  a given  crack  length. 

4.3.1  Average  Strain  Energy  Release  Rate  Using  Area  Method 

From  the  moire  patterns,  we  can  accurately  document  the  strain 
distributions  and  relative  displacements  around  the  crack  tip  as  well  as  the  exact 
location  of  the  crack  tip  which  is  the  key  to  determine  the  crack  length  accurately. 
A simple  method  to  obtain  the  average  Mode  II  fracture  toughness  is  to  calculate 
the  area  under  the  load  displacement  diagram  and  then  use  the  following 
equation 


G 


n 


A U 
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Where  AU  is  the  total  energy  for  crack  propagation  determined  from  the  load- 
deflection  diagram  and  AA  is  the  crack  area  created. 

Two  diagrams  of  the  load  vs.  displacement  of  the  three  point  bending 
setup  are  shown  in  Figure  4-4.  One  can  see  when  the  load  reached  about  489  N 
(110  pounds),  the  slopes  of  the  curves  started  to  decrease  and  the  crack  began 
to  propagate.  Due  to  crack  propagation,  the  flexure  modulus  of  the  specimens 
gradually  reduces.  After  the  load  reached  about  300  pounds,  the  specimens 
were  unloaded  and  other  unloading  curves  were  obtained.  The  area  between 
loading  and  unloading  curves  represents  the  total  energy  for  crack  propagation. 
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Figure  4-5  shows  the  load  vs.  displacement  diagram  by  using  the 
cantilever  beam  setup.  The  displacement  diagrams  of  these  two  specimens 
match  well.  When  the  load  reached  about  40  pounds,  the  crack  started  to 
propagate  for  both  specimens  and  the  slopes  of  the  curves  began  to  gradually 
decrease. 

This  cantilever  beam  setup  gave  us  much  longer  final  crack  lengths  than 
that  of  the  three-point  bending  setup  and  the  maximum  crack  lengths  of  the 
specimens  vary  from  1.75  to  2.75  inches.  It  is  possible  to  obtain  the  fully 
developed  bridging  zone  by  using  the  cantilever  beam  setup  and  the  maximum 
value  of  Gn.effect  can  be  obtained  accordingly. 

Figure  4-6  gives  the  comparison  of  average  stain  energy  release  rates 
from  the  two  different  setups.  The  Gn  of  specimens  1 and  2 using  the  three  point 
bending  setup  is  little  bit  lower  than  the  Gn  of  specimens  3,  4 and  5 using  the 
cantilever  beam  setup.  The  reason  is  that  the  maximum  crack  lengths  of 
specimens  from  the  three  point  bending  setup  are  much  shorter  than  the  ones 
from  the  cantilever  beam  setup.  The  more  stitches  that  are  involved,  the  more 
energy  is  needed  to  cause  crack  propagation.  The  Gn-ave  increases  about  3 to  4 
times  due  to  stitching  compared  to  the  Gn-aVe  of  an  unstitch  specimen, 

4.3.2  Determining  the  Stress  Intensity  Factor  by  Using  High  Sensitivity  Moire 

Interferometry  Technique 

Using  the  area  method,  one  can  only  obtain  the  “average  Ghc”.  In  order  to 
calculate  the  instantaneous  GhC,  one  needs  to  measure  the  crack  length  as  well 
as  the  load  precisely.  Due  to  stable  crack  propagation,  the  load  can  be  read 
from  the  computer  at  any  instant.  Then  we  hold  this  load  and  take  pictures  of  the 
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fringe  patterns.  From  the  moire  pattern,  the  exact  location  of  the  crack  tip  can  be 
determined  and  hence  the  crack  length  can  be  measured  precisely. 

4.3.2. 1 Stress  intensity  factor  for  orthotropic  brittle  materials 

For  orthotropic  materials  under  plane  stress  or  plane  strain  conditions,  the 
relationship  between  stresses  and  strains  can  be  expressed  as 


^x=aucrx+an(7y+a,(lTxy 
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The  equilibrium  equations  for  plane  stress  and  plane  strain  are 
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The  equilibrium  equations  will  be  satisfied  if  the  stress  function  U(x,y)  is 
expressed  as: 
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Substituting  for  equation  (4)  into  the  compatibility  equation 
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The  governing  differential  equation  of  plane  stress  and  plane  strain  of 
orthotropic  materials  can  be  obtained  [29] 
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Defining  the  operators  Oy  (j  = 1,2,  3,  4)  as 

0=1.2,  3, 4)  (7) 

oy  ox 


The  governing  equation  in  U(x,y)  becomes 

D1D2D3D4U(x,y)=0  (8) 


The  jUj  are  the  roots  of  the  characteristic  equation 

«n Hj  ~ 2 aX6Hj  + (2 an  + ««)///  - 2 a26//y  +a22=0  (9) 

The  roots  are  either  complex  or  purely  imaginary  and  cannot  be  real  and 
can  be  expressed  as 


= a,  + //?, , s2  — cc2  + ifi2 , 

s2  =ax-  ij3x , s4  =a2-  ij32 
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The  stress  function  U(x,y)  can  be  expressed  in  the  form 

U(x,y)  = 2Re[U1(z1)+U2(z2)]  (11) 

Let  new  functions  <f>{zx)  = dU/dz-,,  y/{z2)  = dU2/dz2 
The  stress  and  displacements  components  can  be  expressed  as  [29] 


ax  = 2Re[sI2^'(zi)  + s2  V(z2)] 
av=2R  e[^(z,)+V(z2)] 
rtv  =-2Re[5i^(Zi)  + 52^'(z2)] 


ux  =2Re[(a1151  +an  -a16J, )</>(zx)  + {ans2  +al2  - a16s2)y(z2)] 

XZ2)] 


u = 2 ReK-12*9'2  + ^ )(j){zx ) 4-  2 -t  a21--a-^sl 


(12) 


For  pure  Mode  II 
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The  relationship  between  Gn  and  KM  for  the  orthotropic  materials  is: 


r*  V 2 \ta22\M2  , 2a)2  + n66 

gi/=kj/-^J( ) + — 7 

\2  \|  rijj  2(3|| 


(14) 


4. 3. 2. 2 Results  from  moire  fringe  patterns  of  both  setups 

Figure  4-7  shows  a series  of  moire  patterns  under  various  loading 
conditions.  The  closer  the  distance  to  the  crack  tip,  the  denser  the  fringe  pattern. 
These  fringe  patterns  provide  rich  information.  For  instance,  the  location  of  the 
crack  tip  can  be  determined  accurately  and  we  can  obtain  the  strain  distribution 
around  the  crack  tip  and  the  relative  displacement  between  any  two  points.  One 
can  see  when  the  load  is  79.9  lbs,  the  crack  tip  is  behind  the  reference  line  and 
the  crack  is  still  in  the  unstitched  zone.  As  the  load  reaches  138.9  lbs,  the  crack 
tip  is  ahead  of  the  reference  line  and  the  crack  propagates  into  the  stitched  zone. 
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When  the  load  reaches  189.0  lbs,  the  crack  tip  is  far  ahead  of  the  reference,  the 
crack  passes  through  numerous  stitches. 

Fringe  patterns  from  the  cantilever  beam  setup  are  shown  in  Figure  4-8. 
The  stitched  zone  is  ahead  of  the  reference  line  and  unstitched  zone  is  behind 
the  line.  Similar  to  Figure  4-8,  the  crack  tip  is  behind  the  reference  line  when  the 
load  is  31.6  lbs  and  crack  tip  is  right  on  the  line  as  the  load  reaches  31.6  lbs  -- 
that  means  31.6  lbs  is  the  initial  critical  load  for  this  specimen. 

The  crack  surface  displacements  behind  the  crack  tip  of  the  three  point 
bending  setup  are  shown  in  Figures  4-9  - 4-11.  From  Figure  4-9,  it  may  be 
noted  that  the  absolute  values  of  the  upper  and  lower  surface  displacements  are 
equal  for  smaller  loads  (79.9  lbs),  but  they  differ  significantly  as  the  loads  reach 
larger  values  (138.9  lbs  and  189  lbs).  This  may  be  due  to  the  asymmetry  in  the 
specimen  about  the  delamination  plane  or  the  stitching  effect,  which  needs 
further  investigation.  Since  the  density  of  fringe  pattern  is  extremely  high  in  the 
region  very  close  to  the  crack  tip,  it’s  very  difficult  to  obtain  the  displacement  in 
this  area.  Fortunately,  when  one  calculates  the  stress  intensity  factor,  it’s  not 
necessary  to  select  a position  inside  0.5  mm  (0.02  inch)  of  the  crack  zone. 
Similar  results  can  be  obtained  from  the  cantilever  beam  setup. 

4,4  Stress  Intensity  Factor  Kn-namntand  Kn_pffand  Strain  Energy  Release  Rate 

G||-narent  and  G||_fiff 

Usually,  only  the  value  of  Kn  or  GM  is  enough  to  characterize  the  mode  II 
fracture  toughness  of  a material.  However,  for  stitched  composites,  as  a crack 
propagates,  stitches  in  the  bridging  zone  still  partially  connect  to  the  two  crack 
surfaces.  The  value  of  KM  or  Gn  is  not  constant  and  is  directly  related  to  the 
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properties  of  stitches  in  the  bridging  zone.  Here,  author  proposes  two 
parameters,  namely,  Gn-parent  and  Gn.eff  (or  Kn.Parent  and  KM_eff)  to  characterize  the 
fracture  toughness  of  stitched  composites.  Gn-parent  is  one  of  the  intrinsic 
properties  of  the  parent  material  and  should  be  independent  of  the  properties  of 
stitches.  However,  Gn_eff  is  not  only  related  to  the  parent  material  but  also 
depends  on  unbroken  stitches  in  the  bridging  zone.  The  comparison  of  Gn-parent 
and  Gn-eff  will  show  the  effectiveness  of  stitches  in  laminated  composites. 

For  the  purpose  of  computing  KM  using  Eq.  (13)  we  need  only  the  relative 
displacements  between  the  top  and  bottom  crack  surfaces.  The  relative 
displacements  are  compared  for  two  different  loads  (139  and  189  lbs)  in  Figure 
4-12.  It  appears  that  the  relative  displacements  vary  significantly  with  the  loads. 
However,  as  one  moves  closer  to  the  crack  tip,  say  less  than  1 mm  (0.04  inch), 
the  relative  displacement  vs.  distance  curves  for  two  loads  coalesce.  Since  we 
are  using  linear  elastic  fracture  mechanics  (LEFM)  to  predict  the  stress  intensity 
factor  Kn-parent,  one  should  pick  the  position  which  is  as  close  to  the  crack  tip  as 
possible.  The  actual  locations  selected  to  calculate  Kn.Parent  are  away  from  the 
process  zone  but  still  within  the  K-zone.  Further,  the  coalescing  of  the  relative 
displacements  for  those  two  different  loads  indicates  that  the  stress  intensity 
factor  Kn  is  fairly  constant  for  these  loads  which  implies  that  the  crack 
propagation  occurs  at  more  or  less  the  same  Kn-Parent  for  the  crack  lengths 
considered.  It  should  be  mentioned  that  there  is  a relationship  between  KM  and 
Gn  as  shown  in  equation  (14).  We  can  either  use  Kn  or  Gn  to  represent  the  mode 
II  fracture  toughness  of  a material. 
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According  LEFM  the  strains  at  the  crack  tip  should  be  singular,  but  the 
numerical  differentiation  of  displacements  yields  finite  strains  at  the  crack  tip. 
These  values  depend  on  the  yield  stress  of  the  material,  the  finite  curvature  of 
the  crack  tip  and  also  the  location  of  adjacent  stitches  with  respect  to  the  crack 
tip. 

Following  are  the  material  properties  of  the  specimens  used  in  both 
setups. 

Ei=15.06x \06  psi,  E2=E3=1 .60x\06  psi,  v12=0.34,  Gi2=0.80x\0b  psi . 

The  layers  are  orientated  at  45  degrees.  For  the  plane  stress  the  roots  of 
the  characteristic  equation  (9)  can  be  obtained  as 

Si  = 0.3033  - 0.9529i,  S2  = -0.8924  - 0.451 3i; 

S3  = 0. 3033  + 0. 9529i,  S4  = -0. 8924  + 0.451 3i. 

For  the  plan  strain,  the  roots  are 

Si  = 0.3541  - 0. 9352i,  S2  = -0.8952  - 0.4457i; 

S3  = 0. 354 1 + 0. 3529i,  S4  = -0. 8952  + 0. 4457i. 

The  instantaneous  KnC  and  Ghc  can  be  calculated  by  equations  (13)  & 
(14).  Both  are  assumed  in  the  plane  strain  condition.  Figure  4-13  shows  strain 
energy  release  rate  vs.  crack  length  for  laminated  composite  with  stitching 
density  of  40  stitches  per  square  inch.  GM  is  a function  of  the  positions  behind 
the  crack  tip  and  will  reach  a constant  value  as  the  positions  are  close  enough  to 

the  crack  tip.  This  value  is  called  Gn-Parent  and  is  around  3 or  4 (Ib/in).  When 

P=138.9  lbs,  only  one  row  of  stitches  is  in  the  bridging  zone,  the  value  of  Gn.eff  is 
around  10.8  Ib/in  as  the  distance  behind  the  crack  tip  is  more  than  2.5  mm  (0.1 
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inch).  However,  when  P = 189.0  lbs,  there  are  more  than  two  stitches  in  the 
bridging  zone  and  Gn.eff  reaches  a constant  value  18.0  Ib/in  as  the  distance 
behind  the  crack  tip  is  more  than  3.8  mm  (0.15  inch).  This  indicates  that  the 
more  stitches  that  are  present  in  the  bridging  zone,  the  bigger  the  value  of  Gn.eff. 
As  the  bridging  zone  is  “fully  developed”,  which  means  the  bridging  zone  reaches 
a constant  value,  it  will  only  shift  the  location  as  the  crack  continues  to  propagate 
and  Gn-eff  reaches  its  maximum  value.  Figure  4-14  provides  the  strain  energy 
release  rate  vs.  crack  length  for  the  specimens  with  a stitching  density  of  20 
stitches  per  square  inch.  The  results  from  Figure  4-14  are  from  the  cantilever 
beam  setup.  Using  this  setup,  a much  longer  crack  length  can  be  achieved.  The 
crack  length  corresponding  to  P=62.8  lbs  is  about  one  inch  and  at  this  load,  the 
bridging  zone  is  more  or  less  fully  developed.  One  can  tell  from  Figure  4-14  that 
the  maximum  value  of  Gn.eff  is  about  14  (Ib/in)  while  the  value  of  Gn.Parent  is  about 
3.5  (Ib/in).  Gh  increases  about  4 times  for  laminated  composite  with  stitching 
density  of  20  stitches  per  square  inch. 

A previous  study  by  Sharma  and  Sankar  [6]  found  that  for  unstitched 
graphite/epoxy  laminates  Gnc  was  600  N/m  (3.42  Ib/in).  Thus  it  is  clearly  evident 
that  stitching  increases  the  fracture  toughness  by  factors  of  four  to  six  in  the 
specimens  with  stitching  density  of  20  and  40  stitches  per  square  inch 
respectively. 

4.5  Conclusions 

Moire  interferometry  was  used  to  determine  the  displacement  field  in  a 
stitched  graphite/epoxy  ENF  specimen.  From  the  moire  fringe  patterns  the 
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relative  crack  surface  displacements  behind  the  crack  tip  were  determined.  The 
Mode  II  stress  intensity  factor  was  estimated  using  the  linear  elastic  fracture 
mechanics  theory  for  orthotropic  materials.  Another  estimate  of  fracture 
toughness  was  obtained  using  the  area  method,  which  provides  only  average 
value  over  a certain  crack  length.  Results  indicate  that  high  sensitive  moire 
interferometry  is  very  helpful  in  locating  the  crack  tip  during  the  ENF  test  and  the 
relative  crack  surface  displacement  can  be  used  to  compute  the  instantaneous 
values  of  the  Mode  II  fracture  toughness.  Stitching  is  found  to  increase  Gn-effect 
by  factors  of  four  to  six  for  stitching  density  of  20  and  40  stitches  per  square  inch. 
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Figure  4-1 . Three  point  bending  for  end  notched  flexure  test  setup 
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Moire  interferometry  fringe  patterns  on  the  edge  of  an  ENF  specimen 
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Figure  4-4.  Load  vs.  displacement  diagram  of  the  three  point  bending  setup 
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Figure  4-5.  Load  vs.  displacement  diagram  of  the  cantilever  beam  setup 
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Figure  4-6.  Comparison  of  average  strain  energy  release  rates  (Gn)  -specimenl , 
2 tested  by  using  three  point  bending  setup,  specimen  3,  4,  5 by 
cantilever  beam  setup  - Stitching  density  = 20  stitches  /in2 
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Moire  patterns:  U-field 
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Figure  4-7.  Crack  lengths  corresponding  to  the  different  loads  of  the  three-point 
bending  setup 
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Figure  4-8.  Moire  fringe  patterns  under  various  loading  conditions  of  the 
cantilever  beam  setup 
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Figure  4-9.  Relative  displacements  vs.  positions  of  the  three  point  bending  setup 
when  the  load  = 79.9  lbs 
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Figure  4-10.  Relative  displacements  vs.  positions  of  the  three  point  bending 
setup  when  the  load  = 138.9  lbs 
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Figure  4-1 1 . Relative  displacements  vs.  positions  of  the  three  point  bending 
setup  when  the  load  = 189.0  lbs 
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Figure  4-12.  Comparison  between  two  different  loads  of  the  three  point  beam 
setup 
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Figure  4-13.  Strain  energy  release  rate  vs.  crack  length  - stitching  density  = 40 
stitches  per  square  inch. 


86 


Figure  4-14.  Strain  energy  release  rate  vs.  crack  length  - stitching  density  = 20 
stitches  per  square  inch. 


CHAPTER  5 

NEW  MIXED  MODE  TEST  SETUPS 

Until  now,  there  has  been  no  method  to  test  the  fracture  toughness  in 
mixed  Mode  loading  for  stitched  composites.  Methods  that  have  been  developed 
for  non-stitched  composites  cannot  be  extended  to  test  stitched  composites, 
since  stitching  increases  the  fracture  toughness  beyond  that  which  is  possible  to 
test  using  the  existing  methods.  One  needs  to  apply  a high  bending  moment 
which  causes  the  specimen  to  fail  before  crack  propagation.  Since  mixed  mode 
interlaminar  fracture  toughness  is  one  of  the  most  important  parameters  to 
evaluate,  there  is  a great  need  to  develop  test  methods  to  quantify  this  material 
property. 

The  author  has  developed  three  new  test  methods  that  can  be  used  to  test 
mixed  Mode  fracture  toughness  of  stitched  composites.  As  with  the  Mode  I tests, 
the  concept  utilizes  the  idea  of  neutralizing  the  compressive  stress  in  the 
specimen  therefore  preventing  fiber  buckling.  The  ratio  of  Mode  I and  Mode  II 
can  be  adjusted  to  fully  document  the  fracture  toughness  for  any  conceivable 
loading  condition.  Using  this  setup,  one  can  test  numerous  ratios  of  Mode  I and 
Mode  II  and  plot  the  fracture  toughness  with  respect  to  different  ratios  for  stitched 
composites.  These  values  can  provide  a practical  way  to  model  or  design 
stitched  composites  under  different  loading  conditions.  The  test  method  will  be 
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introduced,  load  vs.  displacement  plots  will  be  given  and  the  fracture  toughness 
of  numerous  ratios  will  be  presented. 

5.1  Introduction 

The  author  has  developed  pure  Mode  I and  pure  Mode  II  test  methods  for 
stitched  composites.  However,  in  practical  applications,  it’s  rare  to  encounter 
pure  mode  I or  pure  mode  II  loading  conditions  since  it  is  typical  to  have  a 
combination  of  Mode  I and  Mode  II  loadings.  Until  now,  no  standard  mixed  Mode 
test  was  available  for  laminated  stitched  composites.  Some  researchers 
suggested  to  use  the  MMF  specimen  (Figure  5-1)  [34]  to  test  mixed  Mode 
fracture  toughness  for  regular  laminated  composites.  In  this  setup,  the  ratio  of 
Mode  I and  Mode  is  fixed  at  1.33  according  to  linear  beam  analysis.  Other 
researchers  proposed  the  modified  mixed  Mode  bending  setup  (Figure  5-2)  [35- 
37]  for  different  ratios  of  Mode  I and  Mode  II.  For  this  setup,  numerous  ratios  of 
Mode  I and  Mode  II  can  be  achieved  by  simply  moving  P to  different  positions. 
These  two  setups  are  very  effective  for  unstitched  composites.  Unfortunately, 
both  methods  cannot  be  extended  to  stitched  composites  — especially  heavily 
stitched  composites.  This  is  due  to  the  larger  bending  moment  required  to 
initiate  the  crack  and  cause  crack  propagation  for  stitched  laminated  composites. 
Typically,  one  arm  of  a specimen  buckles  on  the  compressive  side  under  the 
high  bending  moment  before  the  crack  initiates. 

5.2  Test  Procedures 

According  to  our  analysis,  a much  higher  bending  moment  is  required  to 
initiate  the  crack  and  cause  crack  propagation  for  stitched  composites.  In  order 
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to  avoid  buckling,  one  needs  to  substantially  reduce  the  maximum  compressive 
stress.  To  maintain  the  high  bending  moment  while  avoiding  buckling,  one  can 
apply  additional  tensile  force  to  neutralize  the  maximum  compressive  stress.  For 
Mode  I testing,  we  developed  a robust  fixture  which  can  apply  both  transverse 
and  tensile  forces.  The  transverse  force  is  for  opening  the  specimen  and  the 
tensile  force  is  for  reducing  the  maximum  compressive  stress.  Using  this  fixture, 
we  successfully  tested  high-density  stitched  composites  [22-23]  and  measured 
the  Mode  I fracture  toughness.  Compared  to  the  Mode  I test,  the  two  transverse 
forces  in  the  mixed  Mode  test  are  not  equal.  After  some  iteration,  we  developed 
three  new  test  setups  which  can  be  used  to  test  numerous  ratios  of  G|  and  GM 
5.2.1  Novel  Setup  I 

The  key  elements  of  this  fixture  are  the  grips.  Figure  5-3  shows  the 
schematic  drawing  of  the  grip  and  the  gripping  area  of  a specimen.  There  is  a 
Telflon  precrack  at  the  end  of  a specimen  to  simulate  the  initial  crack.  A picture 
of  the  assembly  is  shown  in  Figure  3-4  (Chapter  3)  and  the  profile  of  the  grips 
matches  the  shape  of  the  notches  of  a specimen.  The  grips  can  be  used  for 
specimens  ranging  in  thickness  from  0.125  inch  to  0.300  inch  and  the  grips  can 
hold  both  transverse  and  vertical  force. 

For  this  setup  shown  in  Figure  5-4,  if  a specimen  is  oriented  in  the  vertical 
position  (a=p),  this  setup  is  for  pure  Mode  I since  the  two  transverse  components 
(Fix  and  F2x)  are  equal.  But  as  the  specimen  rotates  as  shown  in  Figure  5-4,  the 
ratio  of  F1x  and  F2x  changes  as  a function  of  a and  p.  By  rotating  the  specimen, 
one  can  obtain  numerous  ratios  of  G|  and  GM.  As  mentioned  before,  the  tensile 
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force  neutralizes  compression.  There  is  a limitation  for  the  maximum  angle  a 
specimen  can  rotate.  If  the  angle  is  beyond  this  limit,  one  arm  of  the  specimen 
will  fail  before  crack  initiation.  In  this  case,  the  tensile  force  is  not  adequate  to 
reduce  the  maximum  compressive  stress  to  avoid  buckling.  For  our  specimens, 
the  maximum  angle  is  around  15  degrees. 

Figure  5-5  shows  photographs  of  the  test  performed  at  a loading  angle  at 
14°.  Using  this  setup,  the  ratio  of  Mode  I and  Mode  II  can  be  adjusted  from 
infinity  to  7.  It  should  be  mentioned  that  due  to  the  deformation  of  the  specimen, 
the  ratio  of  Mode  I and  Mode  II  will  change  slightly  during  a test.  Additionally,  the 
crack  stability  depends  on  the  ratio  of  Mode  I and  Mode  II  and  the  load  vs. 
displacement  diagram  (Figure  5-6)  indicates  the  crack  is  unstable  in  this  test  - 
the  total  number  of  peak  loads  is  less  than  that  of  the  rows  of  stitches. 

5.2.2  Novel  Setup  II 

Novel  setup  II  can  be  achieved  by  only  gripping  one  arm  of  a specimen. 
This  loading  condition  can  be  easily  decoupled  into  the  pure  Mode  I and  pure 
Mode  II  (Figure  5-7).  The  ratio  of  Mode  I and  Mode  II  is  fixed  - no  matter  how 
the  specimen  is  positioned.  There  are  two  major  deficiencies  of  this  setup  - the 
crack  is  unstable  and  the  ratio  of  Mode  I and  Mode  II  is  fixed.  Figure  5-9  shows 
as  the  load  reaches  the  first  peak  load,  significant  energy  releases  and  many 
stitches  break  simultaneously.  If  one  wants  to  obtain  the  instantaneous  Gmax 
(mixed  Mode  critical  strain  energy  release  rate),  this  setup  is  not  a good  choice 
because  it’s  very  difficult  to  capture  the  instantaneous  load  with  respect  to  an 
exact  crack  length. 
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The  ratio  of  Mode  I and  Mode  II  of  this  setup  is  about  1.33  and  neither 
Mode  I nor  Mode  II  is  dominant.  After  observing  broken  stitches  (Figure  5-8),  the 
author  found  they  were  roughly  oriented  at  45°  with  respect  to  the  two  crack 
surfaces  which  indicated  that  the  combination  of  Mode  I and  Mode  II  caused  the 
crack  to  propagate. 

5.2.3  Novel  Setup  III 

On  the  loading  fixture,  the  bar  has  a series  of  holes  (Figure  5-10)  which 
correspond  to  numerous  ratios  of  mixed  Mode.  At  the  right  end  of  the  bar  a 
bearing  reacts  the  horizontal  force  which  balances  the  loading  fixture.  The  free- 
body-diagram  in  Figure  5-10  shows  that  the  ratio  of  Mode  I and  Mode  II  depends 
on  the  loading  positions.  By  changing  the  loading  positions,  one  can  obtain 
different  ratios  of  Mode  I and  Mode  II.  The  ratio  of  Mode  I and  Mode  II  can  vary 
from  infinity  to  zero.  The  fixture  can  be  self-balanced  by  adding  a weight  on  the 
left  side  of  the  bar  to  account  for  errors  caused  by  the  unbalanced  gravity  of  the 
fixture. 

Figure  5-1 1 shows  a series  of  photographs  which  document  the  fracture  of 
a stitched  specimen.  Due  to  mixed  Mode  loading,  the  specimen  bends  to  the  left 
while  at  the  same  time  started  to  open  up.  The  effect  of  Mode  II  caused  the 
specimen  to  bend  to  the  left  and  the  effect  of  Mode  I opened  the  two  arms  of  the 
specimen.  The  combination  of  Mode  I and  Mode  II  causes  crack  propagation. 
Figure  5-12  shows  that  the  total  number  of  rows  of  stitches  matches  the  total 
number  of  peak  loads  and  the  crack  propagates  in  a stable  manner. 
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5.3  Results  and  Discussion 

The  load  vs.  displacement  diagrams  of  the  three  setups  are  shown  on 
Figures  5-6,  5-9  and  5-12.  The  load  increases  until  the  first  row  of  stitches  break, 
at  which  time  the  load  drops  suddenly.  Subsequent  rows  of  stitches  are  loaded, 
then  break,  and  the  cycle  continues  as  a saw  tooth.  After  the  crack  length  has 
reached  the  intended  value,  the  specimen  is  unloaded.  The  total  area  between 
the  loading  and  unloading  curves  is  the  energy  for  crack  propagation.  One  can 
use  the  area  method  to  calculate  the  average  strain  energy  release  rate.  For  the 
third  setup,  the  total  number  of  peak  loads  matches  the  total  rows  of  broke 
stitches  indicating  the  crack  is  stable  and  one  can  use  this  setup  to  determine  the 
instantaneous  Gmix. 

Table  5.1.  Mixed  Mode  strain  energy  release  rate  of  various  ratio  of  Mode  I and 

Mode  II 


Ratio  (Gi/Gn) 

infinite 

7.15 

5.33 

1.33 

0 

G/m“  (Ibin/irf) 

42.4 

44.0 

34.3 

19.0 

0 

G„mix  (Ibin/irf) 

0 

6.2 

6.4 

14.2 

10.7 

Gmix  (ibin/irf) 

42.4 

50.2 

40.7 

33.2 

10.7 

Table  5-1  gives  the  values  of  Gimix  and  Gnmix  in  different  ratios  of  Mode  I 
and  Mode  II.  From  Figure  5-13,  one  can  see  the  values  of  strain  energy  release 
rate  of  pure  Mode  I and  pure  Mode  II  are  not  the  maximum  value  of  G|mix  and 
Gnm,x.  The  maximum  values  of  Gimix  and  Gnmix  occurred  when  the  ratios  were 
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equal  to  7.15  and  1.33  respectively.  From  the  test  data,  we  know  that  the  mixed 
Mode  fracture  toughness  is  directly  related  to  loading  conditions. 

5.4  Conclusion 

Three  novel  setups  were  developed  to  measure  mixed  Mode  properties  of 
stitched  laminated  composites.  These  versatile  setups  can  be  used  to  test  pure 
Mode  I and  pure  Mode  II  as  well  as  any  ratio  of  mixed  Mode.  The  fixture  can  be 
self-balanced  to  eliminate  the  error  caused  by  the  weight  of  the  fixture.  For  some 
ratios  of  Mode  I and  Mode  II,  the  crack  is  stable  and  hence  it  can  be  used  to 
determine  the  instantaneous  Gmix  of  stitched  composites.  Additionally,  this  fixture 
can  be  used  to  determine  the  mixed  Mode  fracture  toughness  for  composites 
with  translaminar  reinforcement  and  has  the  potential  to  be  a standard  test 
method  for  composites  with  translaminar  reinforcement. 
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High  compressive  stress  to  cause 
buckling  under  high  bending  moment 


Figure  5-1,  Mixed  Mode  setup  - ratio  of  G|  and  Gn  is  fixed. 
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Figure  5-2.  Modified  mixed  Mode  setup  - ratio  of  G|  and  Gn  is  adjustable 
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Grip 


Gripping  area 
of  a specimen 


Telflon  precrack 


Figure  5-3.  The  grip  and  the  gripping  area  of  a specimen 
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Figure  5-4.  Novel  Setup  I - by  rotating  a specimen  to  obtain  different  ratios  of  G| 
and  Gii 


Figure  5-5.  Series  of  pictures  taken  from  setup  I 
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Figure  5-6.  Load  vs.  displacement  of  setup  I 
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F 


Figure  5-7.  Novel  setup  II  - gripping  one  arm  of  a specimen 
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Figure  5-8.  Pictures  taken  from  setup  II 
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Displacement  (inches) 


Figure  5-9.  Load  vs.  displacement  of  setup  II 
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Figure  5-10.  Novel  Setup  III  - by  changing  the  loading  position  to  obtain 
numerous  ratios  of  G|  and  Gn 
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Figure  5-11.  Pictures  taken  during  a test  of  setup  III 
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Displacement  (inch) 


Figure  5-12.  Load  vs.  displacement  diagrams  of  setup  III. 
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Figure  5-13.  Relationship  of  Gimix  and  Gnmix. 


CHAPTER  6 

FINITE  ELEMENT  ANALYSIS 

The  tests  described  in  the  previous  sections  were  simulated  using  a finite 
element  analysis  program.  The  purpose  of  performing  the  finite  element  analysis 
was  to  identify  the  correct  physical  models  that  represent  the  stitch  interactions 
with  the  parent  laminate  and  stitch  failure  mechanisms.  The  effects  of  changing 
diameter  and  density  of  the  stitch  yarn  can  be  inferred  without  manufacturing  and 
testing  the  actual  component. 

The  FE  modeling  was  performed  at  two  levels,  two-dimensional  models 
and  full  three-dimensional  simulation.  The  information  from  experimental  data 
was  used  in  generating  the  models.  Specifically,  one  is  interested  in  strain 
energy  release  rate  around  the  crack  tip  and  stresses  in  the  stitches  in  the 
bridging  zone.  The  J-integral  was  used  to  calculate  Gparent  (strain  energy  release 
rate  of  the  parent  material)  by  using  a contour  which  doesn’t  include  the  stitches 
in  the  bridging  zone.  Geff  (the  effective  energy  release  rate)  was  also  calculated 
by  the  J-integral  and  for  this  purpose  the  contour  included  the  stitches  in  the 
bridging  zone.  Gparent  is  a material  property  of  the  parent  laminates  and  should 
be  independent  of  stitch  properties.  However,  Geff  not  only  depends  on  the 
properties  of  the  parent  laminates,  but  also  is  related  to  the  stitches  in  the 
bridging  zone.  The  comparison  of  Gparent  and  Geff  will  give  us  the  effectiveness  of 
stitching  in  laminated  composites.  The  results  in  conjunction  with  those  from 
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experimental  studies  are  expected  to  shed  light  on  the  micromechanics  of  failure 
initiation  and  propagation  around  the  stitches 

The  current  model  assumes  that  the  stitch  interaction  in  DCB  and  ENF 
specimens  can  be  modeled  by  linear  and  nonlinear  bar  elements.  The  inelastic 
behavior  can  be  allowed  to  include  the  ploughing  action  of  the  stitches.  In  the  FE 
2D  model,  the  laminates  were  modeled  using  2D  plane  strain  and  2D  shell 
elements  with  reduced  integration  to  reduce  the  hourglass  effect.  Results  from 
the  2D  models  were  compared  with  the  full  3D  models. 

6.1  Finite  Element  Analysis  for  Mode  I 

6.1.1  Specimen 

The  specimen  was  made  up  of  4 stacks  where  each  stack  consisted  of  7 
plies  which  were  oriented  at  [45°/-450/00/900/-45°/450].  The  materials  used  in 
each  stack  have  slightly  different  properties  as  shown  in  Table  6-1  and  Table  6-2. 


Table  6-1 . One  stack  of  materials 


NUMBER  OF 
PLIES 

THICKNESS 

(INCH) 

ORIENTATION 

(DEGREE) 

NAME  OF 
MATERIALS 

1 

0.00633 

45 

AS4-3501-45 

2 

0.00633 

-45 

AS4-3501-45 

3 

0.01285 

0 

AS4-3501-00 

4 

0.007018 

90 

AS4-3501-90 

5 

0.01285 

0 

AS4-3501-00 

6 

0.00633 

-45 

AS4-3501-45 

7 

0.00633 

45 

AS4-3501-45 
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Table  6-2.  Material  properties 


NAME  OF  MATERIALS 

Ei(PSI) 

E2(PSI) 

Vi  2 

G12(PSI) 

G13(PSI) 

G23(PSI) 

AS4-3501-45 

15.04x10s 

1.6x10s 

0.34 

0.8x10s 

0.8x10s 

0.52x10s 

AS4-3501-00 

15.30x10s 

1.6x10s 

0.34 

0.8x10s 

0.8x10s 

0.52x10s 

AS4-3501-90 

14.88x10s 

1.6x10s 

0.34 

0.8x10s 

0.8x10s 

0.52x10s 

The  specimen  was  7 inches  long  and  0.71  inch  wide  using  stitch  of  1600 
denier  Kevlar  where  there  were  two  yarns  in  each  stitch.  The  stitch  density  was 
8x1/5”  which  means  there  were  40  stitches  per  square  inch  where  the  pitch  is  1/8 
inch,  the  spacing  between  two  adjacent  stitch  rows  is  1/5  inch. 

6.1.2  Effective  Modulus  Ef 

For  the  parent  laminates,  the  specimen  was  assumed  to  have  both 
geometric  and  material  property  symmetry  about  the  neutral  surface  and  each 
ply  is  linearly  elastic  with  no  shear  coupling.  In  reality,  there  are  +45°  and  -45° 
plies  in  laminates  and  shear  coupling  should  exist. 

The  effective  modulus  can  be  easily  calculated  using  the  laminated  beam 
theory  as 

Ej=yt  (emv-zu) 

h 7=i 

Here  N is  the  total  number  of  plies  and  Z;  is  the  distance  from  the  neutral  surface 
to  the  outside  of  the  jth  ply.  Ex  is  the  modulus  along  the  x direction  and  can  be 
expressed  as 
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E 


X 


sin7  a - cos2  a + — sin2  a 


E. 


a = 45°  ->  Ex  =2.109x106(^/) 
a = 0°  ->£,  = 15.3x1 06 (psi) 
a = 90°  ->EX  =1.6x10 6 (psi) 


Where  Ei  E2  and  Gi2  are  the  major  modulus,  minor  modulus  and  the  shear 
modulus  of  each  ply,  respectively  and  a is  the  angle  between  the  major  material 
axis  and  global  x axis.  Ex  is  very  sensitive  to  the  angle  a and  only  about  14%  of 
Ex  along  0°  layer  when  a is  equal  to  45°.  The  effective  modulus  (Ef)  of  this 
specimen  is  equal  to  6.931x106(psi)  and  will  be  used  in  our  2D  solid  model. 

6.1.3  J-inteqral,  “parent  G”  and  “effective  G" 

J integral  is  defined  as  follows: 


Here,  U0  is  the  strain  energy  density,  n-i  is  x component  of  the  outer  normal 
vector  along  the  path  shown  in  Figure  6-1  assuming  crack  propagation  along  the 
x direction,  U ) is  the  displacement  along  the  / direction  and  7,  is  the  component  of 
the  surface  force  along  the  / direction. 

In  Chapter  3,  a new  versatile  fixture  was  used  to  test  the  stitched 
laminated  composites.  Due  to  the  stitching,  the  crack  propagates  in  a stable 
manner.  Numerical  methods  were  used  to  solve  the  nonlinear  governing 
equations  and  the  instantaneous  strain  energy  release  rates  were  obtained  for 
various  crack  lengths.  In  that  analysis,  the  nonlinear  governing  equation  ignored 
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crack  bridging  effects  caused  by  unbroken  stitches  in  the  wake  of  the  crack  tip. 
Hence  that  approach  can  be  considered  as  a method  for  determining  the 
effective  fracture  toughness  (Geff)  of  stitched  laminates.  In  order  to  more 
thoroughly  understand  the  effect  of  the  stitches,  one  should  take  a close  look  at 
the  bridging  zone  where  the  stitches  partially  connect  matrixes  and  prevent  crack 
propagation.  Three  parameters  (Geff,  Gparent  and  COD)  can  be  used  to 
characterize  the  stitch’s  effect  on  the  bridging  zone. 

The  analytical  solution  has  at  least  two  drawbacks  - one  is  that  it  cannot 
exactly  model  the  properties  of  the  parent  materials  (each  layer  is  made  of  3D 
orthotropic  materials),  second  is  that  it  cannot  give  the  exact  crack  tip  opening 
displacement  (COD)  around  the  crack  tip  (the  behavior  of  COD  around  the  crack 
tip  is  totally  different  from  that  of  the  position  away  from  the  crack  tip). 

The  J-integral  is  typically  path  independent  for  the  general  case. 
However,  due  to  the  presence  of  stitches  in  the  bridging  zone,  the  two  crack 
surfaces  are  not  truly  separated  since  they  are  still  partially  connected  by  the 
unbroken  stitches.  Hence,  for  stitched  laminates,  J-integral  is  path  dependent 
when  the  bridging  zone  exists.  If  the  path  of  contour  is  chosen  to  be  very  close 
to  the  crack  tip,  J-integral  should  be  equal  to  Gparent  which  is  a parent  material 
property  regardless  of  whether  the  specimen  is  stitched  or  not.  On  the  other 
hand,  if  the  path  of  contour  includes  unbroken  stitches  in  the  bridging  zone,  the 
J-integral  should  be  equal  to  Geff  which  should  be  larger  than  Gparent  and  will 
reach  a constant  value  when  the  contour  incorporates  all  the  unbroken  stitches  in 
the  bridging  zone.  Subsequently,  as  the  crack  continues  to  propagate,  the 
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bridging  zone  will  eventually  reach  a constant  length  and  move  with  the  crack  tip. 
The  dashed  line  of  the  contour  in  Fig  6.1  was  used  to  calculate  Gparent  and  the 
solid  line  of  the  path  which  includes  all  unbroken  stitches  in  the  bridging  zone 
was  used  for  calculating  Geff. 

6.1.4  2D  Solid  Modeling 
6. 1.4,1  Procedure  of  modeling 

After  carefully  observing  Mode  I specimens,  the  author  found  that  all 
stitches  were  broken  and  no  stitch  was  pulled  out.  Test  results  also  show  that 
there  was  “no  ploughing”  action  for  the  Mode  I test  and  the  average  breaking 
strength  for  Kevlar  1600  denier  is  around  78  lbs  (346  N)  with  the  failure  strain  of 
about  5 percent.  There  are  two  yarns  in  one  stitch  therefore  the  average  failure 
strength  of  each  stitch  is  156  lbs  (694  N).  The  same  ultimate  strength  and  strain 
can  be  used  for  all  stitches  and  the  stitch  can  be  simply  modeled  as  the  elastic 
material. 

The  parent  laminates  are  made  of  3D  orthotropic  materials  and  the 
material  properties  are  strongly  dependent  on  the  fiber  orientation  and  there  is 
some  shear  coupling  effect  in  the  laminates.  In  this  model,  instead  of  modeling 
the  detailed  behavior  of  each  individual  layer,  “effective  Ef”  is  used  to  represent 
the  flexure  modulus  of  the  structure.  It  is  assumed  that  this  model  can  be  used 
to  accurately  calculate  the  displacement  of  the  laminates  and  the  strain  and 
stress  in  the  stitches.  However,  it  cannot  be  used  to  obtain  the  stress  distribution 
in  each  ply.  One  should  be  cautious  when  attempting  to  use  this  model  to 
calculate  the  J-integral  since  it  is  directly  related  to  the  strain  energy  density  and 
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the  gradient  of  the  displacement.  Later  the  J-integral  acquired  from  this  2D  solid 
model  will  be  verified  by  using  a full  3D  simulation. 

Figure  6-2  shows  the  loading  condition  of  the  Mode  I setup  where  all 
forces  transferred  from  the  loading  fixture  are  moved  to  cross  section  A.  As  the 
crack  propagates,  the  crack  length  L5,  measured  to  the  crack  tip  at  B,  increased 
accordingly.  Due  to  stable  crack  propagation,  the  total  number  of  peak  loads 
matches  the  total  number  of  broken  rows  of  stitches  (Figure  6-3).  In  this  model, 
the  author  used  the  critical  loads  as  well  as  the  corresponding  crack  length  from 
experiment  as  input  data.  Table  6-3  lists  the  10  peak  loads  and  the 
corresponding  10  crack  lengths  from  one  of  the  Mode  I test. 


Table  6-3.  Input  data  from  the  experiment  and  analysis 


Pi 

P3 

Ps 

Ps 

Pio 

P12 

P 14 

Pl6 

Pis 

P20 

Loads  (lbs) 

714 

835 

1017 

1137 

1218 

955 

1054 

1140 

1076 

1140 

Fx  (lbs) 

181.1 

206.3 

244.4 

266.2 

279.6 

217.4 

235.2 

251.3 

232.5 

243.1 

Fy(lbs) 

357 

417.5 

508.5 

568.4 

609.3 

477.6 

526.9 

570.2 

537.8 

570.2 

Mz  (lb-ins) 

118.3 

83.1 

45.8 

6.1 

-21.7 

-21.9 

-49.5 

-69.3 

-76.6 

-89.1 

Crack  length  h5  (in) 

0.61 

0.86 

1.11 

1.45 

1.74 

1.99 

2.24 

2.49 

2.74 

2.99 

It  should  be  mentioned  that  the  peak  loads  (Pi  - P2o)  and  the  crack 
lengths  are  obtained  directly  from  the  test  data.  The  loads  (Fx,  FY,  Mz)  at  the 
cross  section  A)  are  calculated  by  an  iterative  method  based  on  the  geometry  of 
the  loading  fixture  and  the  properties  of  the  specimen. 
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Figure  6-4  shows  the  schematic  drawing  of  the  2D  solid  model  of  the  new 
Mode  I test.  The  mesh  around  the  crack  tip  is  much  denser  than  that  of  other 
areas  and  in  this  case,  more  accurate  COD  (crack  opening  displacement)  can  be 
obtained  around  the  crack  tip.  The  stitches  are  connected  on  the  top  and  bottom 
surfaces  of  the  specimen.  The  position  of  the  first  row  of  stitches  is  determined 
under  the  condition  - that  is  under  one  of  the  peak  loads,  the  first  row  of  the 
stitches  is  intended  to  reach  the  failure  strain  or  stress.  One  can  move  the 
position  of  stitches  back  and  forward  in  order  to  find  the  exact  locations  of  the 
stitches  in  the  bridging  zone.  The  bending  moment  is  applied  through  the 
equivalent  couple  forces  FM  (Mz  = FM  d).  The  corresponding  location  of  the  each 
peak  load  comes  directly  from  the  test  data.  In  the  model,  there  are  a total  20 
peak  loads  which  correspond  to  20  locations  of  the  broken  rows  of  the  stitches. 

6.1 .4.2  Results  from  2D  solid  model 

Figure  6-5  shows  the  principal  strain  distribution  under  the  peak  load  Pi 
(714  lbs)  and  the  crack  length  h5  (0.61  inch).  There  are  two  rows  of  stitches  in 
the  bridging  zone.  The  maximum  strain  of  the  stitches  in  the  bridging  zone  was 
5.34%  which  is  approximately  equal  to  the  failure  strain  of  a stitch. 

The  crack  opening  displacement  (COD)  vs.  position  is  shown  in  Figure  6- 
6.  After  a more  detailed  view  in  the  vicinity  of  the  crack  tip,  one  can  see  that 
CODs  are  almost  the  same  (Figure  6-6)  under  various  loading  conditions.  The 
stress  intensity  factor  is  directly  related  to  COD  which  indicates  that  G|.parent  is 
approximately  constant  under  various  peak  loads.  G|.parent  is  an  intrinsic  property 
of  the  parent  material  and  this  model  verifies  our  assumption.  When  the  location 
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is  close  to  the  crack  tip,  COD  vs.  location  can  be  roughly  expressed  as 
$ dose  = C,  Vx  while  moving  away  from  the  crack  tip,  the  relationship  between 
COD  and  the  position  is  approximated  as  Sawav  =C2x2 . Table  6.4  gives  the  G|. 

parent  and  G|.eff  by  using  J-integral  under  various  loading  conditions  and  one  can 
see  that  the  values  of  Gi-parent  are  approximately  same,  however  Gi_eff  varies 
considerably. 


Table  6-4.  Values  of  G|.parent  and  G|.eff  under  various  loads 


Pi 

P3 

Ps 

Ps 

Pio 

P12 

Pl4 

Pie 

P 18 

P20 

loads  (lbs) 

714 

835 

1017 

1137 

1218 

955 

1054 

1140 

1076 

1140 

Gparent  (lb-in/in  ) 

2.7 

3.2 

4.2 

5 

5.6 

3.6 

4.2 

4.7 

3.3 

3.6 

Geff(lb-in/in2) 

93.1 

101 

120.2 

133 

145.1 

111.1 

121.1 

132.4 

122.5 

130 

6.1.5  2D  Shell  Element  Modeling 
6.1. 5.1  Modeling  procedure 


The  model  described  in  6.1.4  uses  the  effective  Ef  to  represent  the  global 
flexural  modulus  of  the  structure  which  provides  realistic  strain  and  displacement 
values.  However,  the  model  in  6.1.4  cannot  be  used  to  obtain  the  stress  field  of 
each  layer.  Due  to  the  symmetric  geometry  of  the  specimen  and  small  shear 
coupling  effect,  2D  shell  elements  are  an  efficient  choice.  The  2D  shell  model 
used  here  is  not  like  the  traditional  plate  model.  From  Figure  6-8,  one  can  see 
that  instead  of  using  the  X-Z  plane,  one  uses  the  X-Y  plane  as  the  mid-plane  for 
2D  shell  elements.  Hence,  there  are  a total  of  28  mid-planes  with  respect  to  28 
plies  in  this  specimen  and  the  “thickness  of  the  shell  element”  is  equal  to  the 
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“width  of  the  specimen”  (along  the  z direction).  Since  each  ply  is  oriented  in  the 
X-Z  plane,  one  must  calculate  Ex  and  Ey  for  each  layer  as  the  input  properties  for 
2D  shell  elements.  Table  6-5  gives  the  material  properties  for  0°,  45°  and  90° 
layers.  Unlike  the  2D  sold  model,  sixth  degree  of  freedom  of  the  2D  shell  model 
is  active  and  the  bending  moment  can  be  directly  applied  to  the  nodes  of  this 
model.  This  is  not  an  exact  model  for  this  structure  but  yields  useful  information 
about  the  stress  field  in  each  ply. 


Table  6-5.  Material  properties  for  shell  elements 


MATERIAL 

X 

LU 

Ey 

Viz 

Gxy 

Gxz 

Gyz 

AS4-3501-00 

1 5.04x1 06 

1.6x106 

0.34 

0.8x106 

0.8x1 06 

0.52x1 06 

AS4-3501-45 

2.1  IxlO6 

1.6x106 

0.34 

0.8x106 

0.8x106 

0.8x1 06 

AS4-3501-90 

1.6x106 

1.6x10s 

0.34 

0.8x106 

0.52x1 06 

0.8x1 06 

6.1. 5. 2 Results  and  discussion 

The  principal  strain  distribution  is  shown  in  Figure  6-9.  The  peak  load  and 
corresponding  crack  length  are  716  lbs  and  0.61  inch  respectively.  Compared  to 
Figure  6-5  of  the  2D  solid  model  under  the  same  loading  conditions,  the  strain 
field  is  almost  identical.  The  maximum  strains  of  the  2D  solid  model  and  the  2D 
shell  model  are  5.18%  and  5.34%  respectively.  The  relative  difference  of  these 
two  maximum  strains  is  about  3%.  If  one  only  wants  to  obtain  strain  and 
displacement  fields,  the  2D  solid  model  is  an  efficient  choice.  Figure  6-10  gives 
the  Mises  stress  distribution  in  the  2D  shell  model.  Under  this  loading  condition, 
the  force  of  the  first  stitch  in  the  bridging  zone  is  164  lbs  and  this  stitch  is  about  to 
break  or  has  already  broken.  The  maximum  stress  on  the  parent  laminates 
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happens  on  zero  degree  layers  with  a value  165  ksi  which  is  less  than  the 
ultimate  strength  (210  ksi)  of  the  AS4-3501  material.  As  the  stress  in  the  stitch 
reaches  the  failure  stress  under  this  peak  load,  the  maximum  stress  of  the  parent 
laminate  is  still  well  below  the  failure  stress. 

This  model  indicates  that  our  new  fixture  can  effectively  break  stitches  and 
maintain  the  parent  laminate  intact.  However,  this  2D  shell  model  cannot  give  us 
the  accurate  J-integral  around  the  crack  tip.  This  is  because  the  mid-plane  used 
in  this  model  is  different  from  that  of  the  traditional  shell  element.  In  order  to 
obtain  G|.Parent  and  G|.eff,  a fully  three  dimensional  model  is  needed  to  simulate  the 
new  Mode  I test. 

6.1.6  Three  Dimensional  Model 

In  this  model,  the  parent  laminate  was  modeled  using  3D-orthotropic 
materials  and  stitches  are  modeled  as  bar  elements.  The  material  properties  and 
orientation  of  the  plies  are  shown  in  Table  6-1  and  Table  6-2. 

The  total  elements  used  in  this  3D  model  are  16004.  3D  solid  elements 
are  used  to  model  the  3D  orthotropic  parent  material  and  ID  bar  elements  are 
used  for  modeling  the  stitches.  Because  the  specimen  undergoes  large 
deformation,  nonlinear  analysis  is  used  to  account  for  the  nonlinear  geometry 
during  static  loading.  Figure  6-11  shows  the  schematic  drawing  of  the  fully  3D 
solid  model  of  the  new  Mode  I test.  The  thickness  of  the  3D  element  is  equal  to 
that  of  the  each  layer  of  the  specimen  and  there  are  a total  28  elements  through 
the  thickness  of  the  specimen.  There  is  only  one  element  through  the  z direction. 
The  positions  of  the  stitches  in  the  bridging  zone  are  determined  by  the  method 
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used  in  the  2D  solid  model.  The  parallel  stitches  are  applied  on  the  “front 
surface”  and  the  “back  surface”  respectively.  Because  the  crack  front  is  a line, 
the  J-integral  is  calculated  through  the  two  lines  of  the  nodes  - one  represents 
the  nodes  on  the  “front  surface”  (the  nodes  in  front  of  the  crack  tip  on  the  front 
surface)  and  another  is  the  nodes  on  the  “back  surface”  (the  nodes  in  front  of  the 
crack  tip  on  the  back  surface). 

After  the  model  was  created,  one  can  simply  manipulate  the  model  to 
move  the  loading  positions  (every  peak  load  has  an  exact  crack  length  from  the 
test)  to  simulate  various  crack  lengths.  This  is  nearly  an  exact  model  in  that  one 
cannot  exactly  model  the  frictional  forces  between  the  stitches  and  the  parent 
laminate  as  well  as  the  frictional  forces  variation  at  different  load  stages 

This  3D  model  can  also  be  used  to  verify  the  2D  model  discussed  in  the 
previous  section.  Figure  6-12  gives  the  comparison  of  COD  in  both  models 
under  the  same  peak  load  P12.  When  the  position  is  very  close  to  the  crack  tip, 
CODs  from  both  models  are  almost  identical.  Moving  away  from  the  crack  tip, 
the  difference  increases  slightly  and  the  maximum  error  is  still  less  than  8%. 
COD  distributions  under  various  peak  loading  conditions  are  shown  in  Figure  6- 
13  and  CODs  are  almost  the  same  for  different  peak  loads  when  the  position  is 
close  to  the  crack  tip.  The  stress  intensity  factor  is  directly  related  to  COD  and 
this  model  shows  that  the  stress  intensity  factors  (K|.parent)  are  roughly  the  same 
for  the  various  crack  lengths. 
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Table  6-6.  Comparison  of  Gi-parent  and  G|.eff  between  two  models 


LOAD 

PI 

P3 

P5 

P8 

P10 

P12 

P14 

P16 

P18 

P20 

3D  Model 

G|-parent 

(lb-in/in2) 

1.5 

1.8 

2.4 

2.8 

3.2 

2.1 

2.4 

2.7 

2.3 

2.5 

2D  Solid  Model 

G|-parent 

(lb-in/in2) 

2.7 

3.2 

4.2 

5 

5.6 

3.6 

4.2 

4.7 

3.3 

3.6 

Relative  Error 

80.0% 

77.8% 

75.0% 

78.6% 

75.0% 

71 .4% 

75.0% 

74.1% 

43.5% 

44.0% 

3D  Model 

G|-eff 

(lb-in/in2) 

67.7 

71.6 

81.1 

89.8 

96.3 

78.1 

83.5 

90.1 

84.1 

88.6 

2D  Solid  Model 

G|-eff 

(lb-in/in2) 

93.1 

101 

120.2 

133 

145.1 

111.1 

121.1 

132.4 

122.5 

130 

Relative  Error 

37.9% 

41.1% 

48.2% 

48.1% 

50.7% 

42.3% 

45.0% 

46.9% 

43.9% 

45.6% 

The  values  of  G|.parent  and  Gi_eff  for  the  2D  model  are  not  a good 
approximation.  The  relative  error  for  G|.parent  and  G|.eff  can  be  as  high  as  88% 
and  50%,  respectively.  G|  is  directly  related  to  the  strain  energy  density  and  the 
2D  model  over-predicts  the  strain  energy.  Figure  6-14  shows  that  Ex  is  only 
2. 1 1 xl 06  psi  for  the  3D  model  for  45°  layers  while  Ef  is  6.11x106  psi  for  the  2D 
solid  model.  If  the  contour  is  chosen  close  enough  to  the  crack  tip,  the  strain 
energy  of  the  3D  model  along  this  contour  is  much  lower  than  the  strain  energy 
of  the  2D  solid  model  under  the  same  strain  condition. 

6.1.7  Diameter  vs.  Density  of  Stitches 

The  stitch  cross-sectional  area  and  the  density  are  the  two  most  important 
factors  when  stitches  are  chosen  as  reinforcement.  Which  parameter,  diameter 
or  density,  has  more  profound  influence  on  the  Gi_eff?  We  use  our  2D  solid  model 
to  vary  these  two  parameters  and  answer  this  question. 
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First,  the  density  of  stitching  is  doubled  and  the  denier  of  Kevlar  is 
maintained.  In  order  to  obtain  the  effective  energy  release  rate  (Gi_eff),  Gi-parent 
should  be  fixed.  We  change  the  loading  conditions  and  the  distance  between  the 
first  stitch  and  the  crack  tip  shown  in  Figure  6-15  until  the  stress  of  the  first  stitch 
in  the  bridging  zone  reaches  the  ultimate  stress  and  G|.parent  is  equal  to  the  value 
obtained  from  the  original  specimen.  After  the  new  LS2  and  the  new  critical 
loading  condition  are  determined,  G|.eff  is  calculated. 

Second,  the  cross-sectional  area  of  a stitch  is  doubled  and  the  density  of 
stitching  is  maintained.  Following  the  same  procedure  described  in  the  previous 
section,  as  the  new  LS2  and  the  new  critical  load  are  determined,  G|.eff  can  be 
calculated  accordingly. 

Table  6-7  provides  results  for  these  two  simulations  and  it  is  apparent  that 
increasing  the  cross-sectional  area  of  the  stitch  has  a more  profound  influence 
on  Gi_eff.  Gi-eff  increases  56.4%  as  the  cross  section  is  doubled  while  G|.eff  only 
increases  28.1%  with  twice  the  stitch  density. 


Table  6.7.  G|.eff  under  two  simulations 


Original  specimen  with  1600  denier 
Kevlar  and  the  stitching  density  8x1/5 

123.5  (lbin/in2) 

Simulation  of  the  specimen  with  the  twice 

158.2  (lbin/in2) 

Relative  increase 

stitch  density 

28.1% 

Simulation  of  the  specimen  with  the 

193.2  (lbin/in2) 

Relative  increase 

doubled  cross-sectional  area  of  stitches 

56.4% 
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6.2  Finite  Element  Model  for  Mode  II 
6.2.1  Modeling  Procedure 

Experiments  revealed  that  for  Mode  II  tests,  there  is  a ploughing  action 
between  the  stitches  and  the  matrix.  Additionally,  the  stitches  seem  to  undergo 
elastic  to  perfect  plastic  behavior  before  breaking.  These  observations  were 
made  during  tests  with  a novel  shear  test  shown  in  Figure  6-16. 

An  extensometer  was  mounted  on  this  novel  shear  loading  fixture  (Figure 
6-16)  to  document  the  relative  displacement  between  the  crack  surfaces.  When 
the  relative  displacement  reached  0.0032  inch  (0.08  mm),  the  load  on  a stitch 
reached  the  maximum  value  67.25  lbs.  The  failure  relative  displacement  of  a 
stitch  is  about  0.008  inch  (0.2  mm). 

The  behavior  of  a stitch  in  Mode  II  loading  is  total  different  from  that  in 
Mode  I loading.  In  the  Mode  I test,  there  is  no  ploughing  action  between  stitches 
and  the  parent  laminate,  however  in  Mode  II,  the  stitches  plough  through  the 
matrix  which  creates  local  damage  of  the  parent  laminate.  Although  the  stitch 
itself  does  not  undergo  plastic  deformation,  the  ploughing  action  can  still  be 
simulated  as  the  plastic  deformation  of  a stitch  (Figure  6-17). 

The  specimen  used  for  the  Mode  II  test  was  6.5  inches  long  and  1 inch 
wide.  Figure  6-18  shows  the  Mode  II  test  setup  and  the  initial  crack  length  of  the 
specimen  was  2.95  inches.  After  loading,  the  two  crack  surfaces  contact  each 
other  and  the  coefficient  of  friction  between  them  was  assumed  to  be  0.15.  High 
sensitivity  moire  interferometry  [24]  was  used  to  determine  the  crack  lengths  and 
the  relative  displacement  between  the  two  crack  surfaces.  Table  6-8  lists  6 peak 
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loads  with  respect  to  6 crack  lengths.  Gn-parent  can  be  determined  by  J-integral 
with  a contour  which  is  in  the  wake  of  the  vicinity  of  the  crack  and  doesn’t  include 
the  unbroken  stitches  in  the  bridging  zone  while  Gn.eff  can  be  calculated  by  a 
contour  which  includes  all  the  unbroken  stitches  in  the  bridging  zone.  As  the 
crack  propagates,  both  Gn.parent  and  Gn.eff  can  be  obtained  under  various  peak 
loads  using  the  J-integral. 


Table  6-8.  Loads  and  crack  lengths  of  Mode  II  test. 


PI 

P2 

P3 

P4 

P5 

P6 

Load  (lbs) 

31.6 

39.7 

46.7 

52.4 

62.6 

62.8 

Crack  length  (inch) 

2.95 

3.03 

3.1 

3.2 

3.54 

3.87 

Figure  6-19  shows  the  schematic  drawing  of  the  full  3D  solid  model  of  the 
new  Mode  II  test.  Like  the  3D  model  for  Mode  I test,  the  thickness  of  the  3D 
element  is  equal  to  that  of  the  each  layer  of  the  specimen  and  there  is  only  one 
element  through  the  z direction.  The  mechanism  of  the  stitches  in  the  Mode  II 
test  is  totally  different  from  that  of  the  Mode  I test.  Due  to  the  ploughing  effect 
between  the  stitches  and  matrix,  one  cannot  directly  connect  a stitch  to  the  top 
and  the  bottom  surfaces  of  a specimen  (in  this  case,  it  cannot  reflect  the  behavior 
of  the  stitches  in  Mode  II  test).  Instead,  “very  short  bar  elements”  oriented  along 
the  crack  surfaces  are  used  to  simulate  the  interaction  between  the  stitches  and 
the  matrix.  The  length  of  the  bar  element  is  0.05  inch  (0.125  mm).  From  the  new 
loading  fixture  (Figure  6-16),  the  stitches’  behavior  in  the  Mode  II  test  can  be 
idealized  as  an  elastic  and  perfect  material  shown  in  Figure  6-17.  It  should  be 
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emphasized  that  the  number  of  the  stitches  in  the  bridging  zone  depends  on  the 
crack  length  and  the  loading  conditions.  For  instance,  under  the  peak  load  P5 
there  are  3 rows  of  the  stitches  in  the  bridging  zone  while  under  the  peak  load  P6 
there  are  4 rows  of  the  stitches  in  the  bridging  zone.  This  model  can  also  be 
used  to  simulate  the  numerous  conditions,  for  instance,  increasing  the  density  of 
the  stitches,  increasing  the  thickness  of  the  stitches  or  even  changing  the 
material  properties  and  the  orientations  of  the  each  layer.  The  previous  full  3D 
model  for  the  new  Mode  I test  was  used  to  simulate  two  conditions  - doubling  the 
density  of  the  stitching  and  doubling  the  cross  section  of  a stitch.  One  can  repeat 
the  same  process  for  this  Model. 

6.2.2  Results  and  Discussion 

Figure  6-20  shows  CODs  under  three  peak  loads  when  the  position  is  very 
close  to  the  crack  tip.  They  are  almost  the  same.  These  results,  again,  verify 
our  assumption  - that  Gn.Parent  is  a parent  material  property  and  should  be 
constant  throughout  a structure.  Two  parameters,  Gn-parent  and  Gn.eff  can  be  used 
to  characterize  the  strain  energy  release  rate  of  stitched  laminates. 

Table  6.9.  Comparison  of  Gii.parent  and  Gn.eff  under  different  loading  conditions 


LOAD 

PI 

P2 

P3 

P4 

P5 

P6 

Gn-parent  (lb-in/in2) 

1.4 

1.5 

1.6 

1.9 

1.6 

1.2 

Gn-eff  (lb-in/in2) 

1.7 

2.6 

3.7 

4.6 

5 

3.9 

Table  6-9  makes  the  comparison  of  Gn.parent  and  Gn.eff  under  different 
loading  conditions.  The  Effective  Mode  II  fracture  toughness  is  strongly 
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dependent  on  the  crack  length  and  the  number  of  stitches  in  the  bridging  zone. 
Usually,  the  more  stitches  that  are  involved  in  the  bridging  zone,  the  larger  the 
Gn-eff.  We  expect  that  as  the  crack  continues  to  grow,  the  first  row  of  stitches  will 
finally  break  and  at  this  moment  the  bridging  zone  will  reach  its  maximum  length. 
After  that,  Gn_eff  will  reach  the  maximum  value  and  will  remain  approximately 
constant  as  the  crack  continues  to  propagate. 
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Figure  6-1 . Contours  of  J-integral 
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Figure  6-2.  Loading  condition  of  the  Mode  I setup. 
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Figure  6-3.  Load  vs.  Displacement  diagram  of  one  specimen  - twenty  peak 
loads  match  twenty  rows  of  the  broken  stitches 


128 


Fx,  Fy,  Mz=FMd  and  a (crack  length)  are  the  input  data 


Figure  6-4.  Schematic  of  2D  solid  model  for  stitched  composites  (the  mesh  of 
the  actual  model  is  much  denser) 


Jr- 


129 


Figure  6-5.  Principal  strain  distribution  of  2D  Model. 
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Figure  6-6.  Crack  opening  displacement  (COD)  vs.  positions 
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Figure  6-7.  COD  vs.  position  close  to  the  crack  tip 
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One  ply  of  laminates 


Figure  6-8.  Schematic  drawing  of  the  mid-plane  for  2D-shell  elements 
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Figure  6-9.  Principle  strain  distribution  in  the  2D-shell  model 
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Figure  6-10.  Mises  stress  distribution  in  the  2D-shell  model 
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Figure  6-11.  Schematic  drawing  of  the  full  3D  solid  model  for  the  new  Mode  I 
(the  actual  mesh  is  much  denser  - the  thickness  of  the  element  is 
equal  to  the  thickness  of  each  layer) 
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Figure  6-12.  Comparison  of  COD  between  the  two  models 
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Figure  6-13.  Crack  opening  displacement  (COD)  vs.  positions  (3D  model) 
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Figure  6-14.  Strain  energy  density  of  the  2D  and  the  3D  model  around  the 
contour 
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Original  specimen  with  Simulation  of  the  doubled 

stitch  density  8x1/5  density  of  stitches 


Figure  6-1 5.  Simulation  of  the  specimen  with  the  twice  stitch  density 
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Figure  6-1 6.  Novel  shear  test  to  determine  the  behavior  of  stitches  under  shear 
loading  conditions. 
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Figure  6-17.  Model  of  stitch’s  behavior  in  Mode  II 


142 


Figure  6-1 8.  Setup  for  the  mode  II  test. 
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one  end  connected  to  top 
surface,  another  end  connected 
to  bottom  surface. 


Figure  6-19.  Schematic  drawing  of  full  3D  solid  model  for  Mode  II  test  (the  actual 
mesh  is  much  denser  - the  thickness  of  the  element  is  equal  to  the 
thickness  of  each  layer) 
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Figure  6-20.  Crack  opening  displacement  for  Mode  II  under  three  different 
loading  conditions 


CHAPTER  7 

SUMMARY  AND  DISCUSSION 

The  standard  DCB  setup  cannot  be  used  to  test  the  Mode  I fracture 
toughness  of  composites  with  medium  to  high  density  stitches.  Although  fiber 
composites  are  strong  in  tension,  their  compressive  strength  is  limited  by  fiber 
micro-buckling,  which  causes  the  specimens  to  fail  before  crack  propagation.  A 
novel  Mode  I test  was  developed  to  test  the  fracture  toughness  of  stitched 
laminated  composites.  It  can  apply  the  extra  tensile  force  required  to  neutralize 
the  maximum  compressive  stress  on  the  outer  surfaces  of  a stitched  specimen. 
Two  analysis  methods  were  developed  to  compute  the  instantaneous  Mode  I 
fracture  toughness  (Gic).  One  is  for  the  preliminary  loading  fixture  and  G|C  can 
be  calculated  as  AW/(Aa*b).  Another  method  utilizes  an  iterative  method  to 
solve  the  nonlinear  governing  equations  to  obtain  the  bending  moment  in  the 
crack  tip  and  Gic  can  be  obtained  as  M2/(bEI).  The  test  method  is  versatile  and 
easy  to  employ 

Two  setups  were  used  to  obtain  the  Mode  II  fracture  toughness  of  stitched 
laminated  composites.  In  both  setups,  high  sensitivity  moire  technique  was  used 
to  document  the  full  displacement  field  as  well  as  the  exact  position  of  the  crack 
tip.  Both  setups  show  stable  crack  propagation  due  to  stitching.  In  the  Mode  II 
test,  the  specimen  undergoes  large  displacement  and  rotations  making  it  difficult 
to  apply  high  sensitivity  moire  technique.  The  large  carrier  pattern  should  be 
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added  to  capture  the  displacement  field  around  the  crack  tip.  The  moire 
technique  can  be  used  to  obtain  the  relative  displacement  between  two  crack 
surfaces,  which  is  directly  related  to  the  stress  intensity  factor.  When  the  position 
behind  the  crack  tip  is  close  to  the  crack  tip,  less  than  0.01  inch  (0.25  mm),  the 
relative  displacement  can  be  used  to  calculate  the  Kparant  (or  Gparent)-  However, 
as  the  position  exceeds  this  range  between  0.01  inch  (0.25  mm)  and  0.25  inch 
(6.4  mm),  the  relative  displacement  can  be  used  to  obtain  the  Keff  (or  Geff).  One 
should  be  careful  when  attempting  to  use  this  relationship  to  obtain  the  Ketf  (or 
Geff)  since  these  values  vary  with  the  location  one  selects.  If  one  wants  to  obtain 
a fully  developed  bridging  zone,  the  second  setup  is  a better  choice  because  it 
can  provide  a much  longer  crack  length. 

In  Chapter  5,  three  new  setups  are  used  to  test  mixed  Mode  fracture 
toughness.  Deficiencies  of  the  first  setup  are  that  there  was  a limitation  of  the 
initial  angles  a specimen  can  be  oriented  (if  the  angle  is  too  big,  the  specimen  is 
going  to  break)  and  also  the  ratio  of  Mode  I and  Mode  II  is  not  constant  during 
the  test.  It’s  difficult  to  identify  the  ratio  of  Mode  I and  Mode  II  since  it  changes 
during  a test.  For  the  second  setup,  the  crack  is  unstable  and  ratio  is  almost 
fixed  (there  is  a slight  change  of  the  ratio  of  Mode  I and  Mode  II  due  to  the  axial 
forces  although  the  contribution  of  the  axial  forces  to  Gic  is  small).  If  one  wants 
to  obtain  the  instantaneous  Gmix,  the  second  setup  is  not  good  choice.  The  third 
setup  is  the  best  choice  to  test  the  mixed  Mode  fracture  toughness  for 
composites  with  the  translaminar  reinforcement.  It  has  many  advantages,  for 
instance,  the  ratio  is  adjustable  (one  setup  can  be  used  to  produce  the  pure 
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Mode  I and  pure  Mode  II  and  any  conceivable  ratio  of  Mode  I and  Mode  II),  the 
ratio  changes  little  during  a test  (due  to  the  unbalanced  forces  applied  from  the 
fixture,  the  specimen  will  tilt  to  one  side  and  hence  the  ratio  will  change  slightly 
during  a test)  and  the  fixture  can  be  self-balanced  to  eliminate  the  error  caused 
by  the  gravity  of  the  fixture.  Due  to  the  lack  of  specimens,  only  two  ratios  of 
Mode  I and  Mode  were  tested  by  using  the  third  setup.  The  author  hopes  in  the 
future  to  test  at  least  6 to  10  different  ratios  using  this  fixture  and  then  based  on 
the  experimental  data,  create  a fracture  criteria  for  composites  with  the 
translaminar  reinforcement. 

FEA  models  proposed  in  Chapter  6 provide  valuable  information  about 
stitched  composites.  The  potential  of  these  models  is  to  obtain  accurate 
behavior  of  the  stitches  during  the  Mode  I and  Mode  II  tests.  The  behavior  of 
stitches  alone  is  totally  different  from  that  of  stitches  inside  the  parent  materials. 
For  the  Mode  I test,  stitches  are  modeled  as  bar  elements  and  stitches  act  like 
elastic  material.  They  are  connected  to  the  top  and  bottom  of  the  parent 
material.  In  reality,  there  exists  the  interaction  between  the  matrix  and  the 
stitches  and  the  strain  and  stress  should  vary  throughout  the  length  of  a stitch.  In 
order  to  more  thoroughly  understand  the  stitch’s  behavior  inside  the  stitched 
composites,  the  author  suggests  using  the  new  fixture  described  in  Chapter  3 to 
pull  out  one  or  two  stitches  from  the  parent  material.  This  test  could  shed  light  on 
the  micromechanics  of  stitch  failure.  Then  based  on  this  test  data,  more 
accurate  model  of  a stitch  in  Mode  I can  be  established.  In  Chapter  6,  only  two 
simulations  (doubling  the  density  of  the  stitching  and  doubling  the  cross  section 
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of  a stitch)  are  used  to  determine  the  fracture  toughness  of  stitched  laminated 
composites.  Many  parameters  can  be  simulated  by  using  these  models,  for 
instance,  changing  the  orientations  of  each  layer,  thickness  of  the  specimen  and 
even  the  properties  of  the  parent  material  and  stitches. 

The  behavior  of  stitches  in  Mode  II  is  based  on  the  new  shear  test  (Figure 
6-16).  An  extensometer  is  mounted  on  the  specimen  to  document  the  relative 
displacements  of  the  two  crack  surfaces.  The  load  vs.  displacement  diagram 
was  obtained  after  the  two  crack  surfaces  are  separated  (the  stitches  still 
connect  the  two  separated  surfaces).  The  stitch’s  behavior  under  Mode  II 
loading  is  totally  different  from  that  under  the  Mode  I loading.  In  Mode  II  test, 
stitches  produce  “shear  forces”  to  prevent  the  two  surfaces  from  sliding  with  each 
other.  The  stitches  also  cause  local  damage  of  the  matrix.  There  is  the  friction 
between  the  two  contact  surfaces  (they  contact  with  each  other  before  and  after 
crack  propagation).  The  coefficient  of  the  friction  may  vary  under  various  loading 
conditions.  It’s  impossible  to  accurately  model  these  effects.  In  our  model,  a 
“very  short  bar  element”  oriented  along  the  two  crack  surfaces  is  used  to 
simulate  the  behavior  of  a stitch  in  the  Mode  II  test.  This  is  not  an  exact  model, 
however  it  certainly  predicts  the  major  behaviors  of  stitches  in  the  Mode  II  test. 
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